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ABSTRACT 

We model the luminosity-dependent projected two-point correlation function of DEEP2 and SDSS 
galaxies within the Halo Occupation Distribution (HOD) framework. From this we infer the rela- 
tionship between galaxy luminosity and host dark matter halo mass at 2 ^ 1 and at z ~ 0. At 
both epochs, there is a tight correlation between central galaxy luminosity and halo mass, with the 
slope and scatter decreasing for larger halo masses, and the fraction of satellite galaxies decreasing 
at higher luminosity. Central L, galaxies reside in halos a few times more massive at z ~ 1 than 
at z ~ 0. The satellite fraction of galaxies more luminous than L» is ^10% at z 1, compared to 
~20% at z ~ 0. We find little evolution in the relation between the characteristic minimum mass 
of a halo hosting a central galaxy above a luminosity threshold and the mass scale of a halo that on 
average hosts one satellite galaxy above the same luminosity threshold, with the latter being 15-20 
times the former. Combining these HOD results with theoretical predictions of the typical growth of 
halos, we establish an evolutionary connection between the galaxy populations at the two redshifts 
by linking z ~ central galaxies to z ~ 1 central galaxies that reside in their progenitor halos, which 
enables us to study the evolution of galaxies as a function of halo mass. We find that the stellar mass 
growth of galaxies depends on halo mass. On average, the majority of the stellar mass in central 
galaxies residing in z ~ low-mass halos ('--^ 5 x lO^^/i~^M0) is the result of star formation between 
z ~ 1 and z ~ 0, while only a small fraction of the stellar mass in central galaxies of high mass 
halos IQ^^^h'^Mo) is the result of star formation over this period. In addition, the mass scale 
of halos where the star formation efficiency reaches a maximum is found to shift toward lower mass 
with time. With appropriately defined galaxy samples at different redshifts, future work can combine 
HOD modeling of the clustering with the assembly history and dynamical evolution of dark matter 
halos. This can lead to an understanding of the stellar mass growth due to both mergers and star 
formation as a function of host halo mass and provide powerful tests of galaxy formation theories. 
In the appendix, we provide a brief discussion of systematic biases related to the assumption of "one 
galaxy per halo" in estimating the mass scale and number density of host halos from the observed 
clustering strength of galaxies. 

Subject headings: cosmology: observations — galaxies: clustering — galaxies: distances and red- 
shifts — galaxies: evolution — galaxies: halos — galaxies: statistics — large-scale 
structure of universe 



1. INTRODUCTION 

The spatial clustering of galaxies encodes useful in- 
formation about their formation and evolution, which 
remain outstanding problems in modern cosmology and 
astrophysics. Large redshift surveys at different epochs 
have enabled detailed studies of galaxy clustering and 
its evolution. In particular, the dependence of clustering 
on galaxy properties such as luminosity and color pro- 
vides fundamental constraints on galaxy formation theo- 
ries. In this paper, we model the luminosity-dependence 
of the galaxy two-point correlation functions at z '--^ 1 
and z ~ measu red for galaxies in the DEEP2 Galaxy 
Redshift Survey (ICoil et al.ll2006a^ and t he Sloan Digi- 
tal Sky Survey fSDSS: (jZehavi et al.|[2005l ) to investigate 
the evolution of galaxies over the last ^7 billion years. 
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To make full use of the galaxy clustering measurements 
to test galaxy formation theories, one needs to extract 
and characterize the information available in the clus- 
tering data in a convenient and informative form. In 
the cold dark matter (CDM) hierarchical structure for- 
mation paradigm, the evolution of galaxies is coupled to 
that of the dark matter halos, defined as roughly spheri- 
cal, virialized regions with over-density about 200 times 
that of the mean background density. The formation 
and evolution of dark matter halos are dominated by 
gravity and, with the help of improving computational 
power and iV-body simulations, can be calculated ac- 
curately for any specified cosmological model. Galaxy 
clustering data are particularly illuminating if one can 
relate galaxies to dark matter halos in a way that results 
in informative tests of galaxy formation theory. This 
can be achieved within the framework of the halo oc- 
cupation distribution (HOD), which describes the sta- 
tistical relati on between galaxies and dark matter ha- 
los (see e.g.. iJing Mo fc Borndll998l: iMa fc Frvll2000l: 
Peacock & Smith! 120001; lSeliak"2000'; 'Scoccimarro et aU 
2001; Bcrhnd & Weinberg 2002; Cooray & Sheth 2002). 
The HOD characterizes this relation in terms of the prob- 
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ability distribution P{N\AI) that a halo of virial mass M 
contains N galaxies of a given type and the relative spa- 
tial and velocity distribution of galaxies and dark matter 
within halos. 

Contemporary large galaxy surveys, such as the SDSS 
(jYork et all [200Q1 and the Two-Degree Field Galaxy 
Redshift Survey r2dFGRS: IColless et alJ [200ih . enable 
detailed studies of the z ~ galaxy population. HOD 
modeling, or the closely rel ated "conditional luminosity 
funct ion" (CLF) method (lYang. Mo. fc van den BoschI 
|2003| ). has been apphed to interp ret clustering 
data from these local su r veys (e.g., [jTng fc Borned 



1998[ iJing. Borne r. fc Sutol 12002 1: 'van den Bosch et all 
20031: [M aghocchet ti fc Porcianil l2GG3: Zehavi et al.ll2005t 
Yang et al.l 120051) . Recently, much effort has also 



been placed on studying the clustering of high red- 
shift gal^jdeSj_from_2: < 1, such as the DEEP2 sur- 
vey ("Davi s et al.ll2003f) . the COMBO-17 survey (Classi- 
fying Obj ects with M edium Band Observations in 17 fil- 
ters; iWolLiEH^'ZOOl), the NDWFS (N O AO Deep Wide- 
Field Survey; Jannu zi fc Deyl 11993). an d the VVDS 
(VIMOS-VLT Deep Survev: ILc Fevre et al. 200 g), up 



to z 



3-5 in many sur v evs (e.g.. [Adelbereer et al.l 



2005allbl: iLee et all l20oi lOuchi et all 12004 12005 



Allen et all l2005t iHildebrandt et all l2005t iDaddi et al l 
20031 ). HOD/CLF modehng has also been used to ex- 



plain the cluste ring properties of these high r e dshift 
galaxies (e.g., 'Bullock. Wechsler. fc S omervilld 120021: 
Moustakas fc Somcrvillc 2002; Yan, Mad gwick. fc White! 
2003; Zheng 2004; Lee ct al. 2006; Hamana et al. l 120061: 
Coorav 2005, 200 6: Coorav fc Ouchi 2006; Conr ov et all 
20061: 1 White et ani2007l : Is'lake et al.ll2007l) . Galaxy clus- 



tering data at different redshifts provide a goldmine of 
information from which to deepen o ur understanding o f 
galaxy evolution. As an example, 'Wh ite et al.l ()2007D 
find evidence for merging or disruption of red galaxies in 
NDWFS from modeling the evolution of their clustering. 
Recent studies of galaxy luminosity functions and two- 
point correlation functions at different redshifts with the 
CLF approach also reveal some interesting evolutionary 
trends, such as a brightening of the luminosity of cen- 
tral g alaxies at a fi xed halo mass toward higher redshifts 
(e.g.. lCooravll2005l) . 

As galaxies evolve, they grow in both stellar and 
dark matter mass. Their stellar component continuously 
changes as the original stars age and new stars form. 
Galaxies can also change their stellar contents and in- 
crease their mass as a result of mergers with other galax- 
ies. Different feedback or pre- heating mechanisms, such 
as those caused by star formation, active galactic nu- 
clei, or the photo-ionizing ultra-violet background, can 
also impact at different stages of a galaxy's life. Var- 
ious evolutionary processes can change galaxy proper- 
ties such as luminosity, color, mass and morphology. To 
observationally establish a connection between galaxy 
populations at different redshifts is therefore not triv- 
ial. Much of the work to identify progenitors or de- 
scendants of certain galaxy types from galaxy cluster- 
ing has been based on rather simple models with limited 
constraining po wer, such as the object - conserving or the 
merg ing model (jMatarrese et al.lll997l : iMoscardini et al.l 
11998^. By contrast, the HOD framework converts the ob- 
served galaxy clustering data to a relation between galax- 
ies and their host dark matter halos, placing galaxies in 



a cosmological context. Using HOD constraints at differ- 
ent redshifts derived from clustering data thus provides 
a physically motivated way to study galaxy evolution. 

The evolution of the HOD reflects how the relation be- 
tween galaxies and halos changes with time, which can 
be used to test galaxy formation models in a more trans- 
parent way. The evolution of dark mat ter halos can be 
solved with analytic approaches (e .g.. iTavlor fc Babull 
200ll [200l iBc nson ct al. 200 l lOguri fc Led 12004 



van den Bosch et al.l 120051 : (Zentncr et al. 2005") and cos- 
mological simulations (e.g.7ISpringcl_ct__al. 2005). Thus, 
the HOD evolution may potentially provide valuable con- 
straints on the evolution of baryons in galaxies (e.g., 
accretion and consumption of gas, formation of stars) 
and put powerful tests on galaxy formation theories. 
Many recent studies of galaxy evolution have mea- 
sured specific galaxy properties (e.g., the star forma- 
tion rate or the star formation history) as a function 
of the mass of the galaxy, where either the stellar mass 
or an estimate of t he dynamical mass is used (e.g.. 
Heavens et al.l l2004| iJimenez et al.l [20051: lJuneau et al. | 



2005; Bundv ct al. 200' 
2006aiib: Jimenez et al 



Cimatti et all 120 06: Er b et all 
20061 iNoeskeet al.i i2007a 



Zheng et al. I l2007fl . The unique aspect of the approaches 



presented and envisioned in this paper that make use of 
HOD modeling of galaxy clustering at different redshifts 
is that we are able to study the evolution of galaxy prop- 
erties as a function of the host dark matter halo mass, 
which is a fundamental parameter in galaxy formation 
and evolution. 

In this paper we perform HOD modeling of the lu- 
minosity dependence of the galaxy two-point correlation 
function for samples at two epochs, the z ~ 1 DEEP2 
galaxies and the z ~ SDSS galaxies, to study the evolu- 
tion of the HOD and its implication for galaxy formation 
models. The structure of the paper is as follows. In § 2 
we briefly describe the DEEP2 and SDSS galaxy samples 
used here. We define our cosmological model and intro- 
duce the HOD parameterization in § 3. The modeling 
results are presented in § 4. In § 5 we make a detailed 
comparison between the HODs at the two epochs. In 
§ 6 we attempt to establish an evolutionary link between 
the two galaxy populations through the growth of dark 
matter halos. This allows us to infer the evolution of 
galaxies as a function of the host halo mass during the 
last ^7 billion years. We also discuss the star formation 
efficiency at the two epochs. Finally, in § 7 we sum- 
marize our results and discuss how to extend the study 
in this paper to a powerful, comprehensive program to 
constrain galaxy formation and evolution theories. For 
a simple estimate, instead of using a full HOD model- 
ing, some applications of galaxy clustering assume one 
galaxy per halo. In the Appendix we use an HOD model 
to investigate the potential systematic errors introduced 
by this assumption. 

2. GALAXY SAMPLES 

2.1. DEEP2 Samples 

We use th e projected two-po int correlation functions 
measured bv lCoil et al.l (f2006a) for luminosity threshold 
galaxy samples from the DEEP2 Galaxy Redshift Sur- 
vey to constrai n the HOD at z ^ 1. The DEEP2 survey 
uses DEIMOS (|Faber et al.ll2003l ) on the 10-m Keck II 
telescope to survey optical galaxies at z '--^ 1 in a comov- 
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ing volume of approximately 5x10^ h^^Mpc^. DEEP2 
has measured redshifts for > 30, 000 galaxies in the red- 
shift range 0.7 < z < 1.45 to a hmiting magnitude of 
Rab — 24.1 over 3 deg^ of the sky. Technical de t ails o f 
the DEE P2 survey can be found in [Pavis et all (|2003l ) 
and Davis. Gerke. fc NewmanI (|2004f l. 

Here we use the four nearly volume-limited samples of 
I Coil et al] (|2006af ). defined by thresholds in rest-frame 
Johnson B absolute magni tude (Mb)- ivT - correc tions are 
calculated as described in Will mer et al.l ()2006[ ). and no 
corrections are made for luminosity evolution. Each sam- 
ple has a minimum redshift of z = 0.75 and a maximum 
redshift of z = 1.0 — 1.2, depending on luminosity, and 
includes 5000-11,000 galaxies. The three brightest sam- 
ples (Mb < -20.5, Mb < -20.0, and Mb < -19.5) are 
volume-limited for blue galaxies, while the faintest sam- 
ple {Mb < —19.0) is missing some fainter blue galax- 
ies at 0.95 < z < 1.0. The samples are not entirely 
volume-limited for red galaxies, due to the Rab selec- 
tion of the DEEP2 survey (for definitions of blue and 
red galaxies and a discussion of this selection effect see 
I Willmer et aT][2006h . By assuming little evolution in the 
red galaxy fraction at z ~ 1 and using a volume-limited 
(but smaller) control sample at a slightly lower redshift, 
we estimate that, for each sample, about 10% of the full 
population is missed because of the red galaxies being not 
volume-limited. We address the impact of this sample se- 
lection on the galaxy evolution we attempt to establish 
in §E1 

2.2. SDSS Samples 

We compare the DEEP2 clustering and modeling re- 
sults to thos e obtained from the SDSS at lower redshi fts. 
The SDSS (|York et al.l [20001 : iStoughton et~an [200l ) is 
an ongoing project that aims to map nearly one quar- 
ter of the sky in the northern Galactic cap, as well as 
a small portion of the so uthern Galactic ca p, using a 
dedicated 2.5m telescope (iGunn et al.ll2006D. A drift- 
scanning mosaic CCD camera ()Gunn et al.lll998D is used 
to image the sky in five photometric bandpasses to a 
limiting magnitude of r ~ 22.5. Objects are selected 
for spectroscopic follow-u p using specif ic algorithms for 
the main galaxy sa mple (IStrauss et al.1 [2002,1 . luminous 
red galaxy sample (lEisenstein et al.ll2001f ). and quasars 
([Richards et al.l [200 2^1. To a good approximation, the 
main galaxy sample consists of all galaxies with Pet- 
rosian magnitude r < 17.77, with a median redshift of 
~ 0.1. The construction of t he large-scale st r uctur e sam- 
ples is described in detail in iBlanton et al.l ([2OO50 . The 
radial selection function is derived from the sam ple se- 
lection criteria using the i^T-corrections of Blanto n et al.l 
((2003a) and an impro ved version of the evo lving luminos- 
ity function model of IBlanton et al.l ([2003b) . The angu- 
lar completeness is characterized carefully for each sector 
(a unique region of overlapping spectroscopic plates) on 
the sky. 

We use here t he clustering results presented in 
[Zehavi et al.l ([2005[ ) for an SDSS sample containing some 
200, 000 galaxies extending over ~ 2500 deg^ on the sky. 
More specifically, we use the projected two-point cor- 
relation function measurements for volume-limited sam- 
ples defined by thresholds in r-band absolute magnitude 
(Mr). These luminosity thres hold samples span a red- 
shift range of 0.02 < z < 0.22. [Zehavi et all ([2005f) have 



also performed detailed HOD modeling of these cluster- 
ing measurements. We build on these results, although 
in order to facilitate a more direct comparison to the 
DEEP2 modeling presented here, we repeat the SDSS 
HOD modeling with a slightly different and more flexi- 
ble suite of models used in this work. 

The SDSS and the DEEP2 surveys have different sam- 
ple selections, with SDSS galaxies selected in rest-frame 
r-band and DEEP2 galaxies in rest-frame B-band. The 
differences in the selection complicate the comparison of 
the HOD modeling results, a point we address through- 
out the paper. 

3. COSMOLOGY AND HOD MODELS 
3.1. Cosmology and Halo Properties 

We adopt a spatially-flat CDM cosmological model 
with Gaussian initial density fluctuations that is consis- 
tent with the determination from the cosmic microwave 
backgroun d. Type la supernovae, and galaxy c l uster- 
ing fe.g.. [Spergel et al.[ [2003[. [2007t [Riess et al.[ [2004 
[Tegmark et al.[[2004t[Abazaiian et al ]|2OO50. Density pa- 
rameters are assumed to be (l^m, ^a, r2b)=(0.3, 0.7, 
0.047), where flm, ^^A, and fib are density parame- 
ters of the matter (CDM-|-baryons), dark energy, and 
baryons, respectively. The mass fluctuation power spec- 
trum is the primordial power spectrum, with the spec- 
tral index assumed to be ris—l, modifled by the trans- 
fer function, which i s computed using the formula given 
by [Eisenstein fc Hu[ (fl99S) with the effect of baryons 
taken into account. It is normalized such that the rms 
fluctuation of the linear density in spheres of radius 
Sh^^Mpc at z=0 is cr8=0.9. The Hubble constant is 
Ho = lOO/i = 70kms~^Mpc"^ 

We define a dark matter halo as an object with mean 
density of 200 times that of the background universe. 
With the assumed cosmological model, properties of dark 
matter halos are calculated based on numerically tested 
fitting formulas. The halo mass function is cq r npute d 
according to the formula given bv [Jenkins et al.[ ([2001h . 
For the large-scale h a lo bia s factor, we adopt the for- 
mula in [Tinker et al.l ([20 05'), which is a modification of 
that given by [Sheth et a l. (2001) and is accurate for a 
large range of cosmological parameters. The density 
distribution of a dark matter halo of mass M is as- 
sum ed to follow the Navarro-Frenk- White (NF W) pro- 
file (iNavarro. Frenk. fc Whitd[l995l [19961 [l99l charac- 
terized by the concentration parameter c(M). For c(M) , 
we use the relation given bv lBullock et~an (|200lh . mod- 
ified to be consistent with our halo definition. 



c(M) 



Co 



1 



■{M/M, 



(1) 



where cq = 11, /3 = —0.13, and is the nonlinear 
mass at z = (Af* = 7.77 x IO^^/i-^Mq for the adopted 
cosmology) . 

3.2. HOD Parameterization 

It has proven to be useful to parameterize the 
HOD by separating contributions from central and 
satellite galaxies, through studies of subhalo HODs 
in high-resolution, d issipationless A^-body simulations 
([Kravtsov et al.[[2004[ ) and galaxy HODs in semi-analytic 
galaxy formation model (SA) and in smoothed particle 
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hydrodynamics (SPH) simulations (jZheng et al.l 120051 ) . 
For galaxy samples defined by lower luminosity thresh- 
olds, a simple parameterization of the mean occupation 
function has three parameters (e . g..lKravtsov et al.|[200^ 
IZheng et all 120051 IZehavi et all l2005f ) — the mean oc- 
cupation function {Ncen{M)) of central galaxies can be 
represented by a step-like function with a characteris- 
tic minimum halo mass Mmin, and the mean occupation 
function {Nsg^t{M)) of satellite galaxies is approximated 
by a power law with the amplitude and slope as two free 
parameters. This simple three-parameter model can cap- 
ture the basic features in the mean occupation function 
predicted by galaxy formation models. 

However, there are some advantages in choosing a 
parametrization more flexible than the simple one. For 
example, at lower halo mass, the mean occupation func- 
tion of satellites predicted by galaxy formation models 
drops faster than the extrapolati o n of the high mas s 



drops taster tnan tne extrapolati o n ot tne m; 
power law (jKravtsov et al.l T2004I : IZheng et~al 



l2005l) . 

Since, in the low-mass range, satellites have a substan- 
tial contribution to the galaxy two-point correlation func- 
tion, with the simple three-parameter model, the best-fit 
power law slope in the satellite mean occupation func- 
tion tends to reflect the effec tive slope in the low-mass 
range (e.g.. lZehavi et ani2005l ). That is, the overall slope 
is adjusted to flt the amplitude of (Ngat) in the low-mass 
range, and it may become a poor description for the slope 
at high halo mass. 

In this paper we adopt a slightly more flexible parame- 
terization with five parameters, motivated by the results 
presented in|Zhcng ct al. (2005). The mean occupation 
function of the central galaxies is a step-like function 
with a soft cutoff profile to account for the scatter be- 
tween galaxy luminosity and host halo mass. That of 
the satellite galaxies is a power law modified by a low- 
mass cutoff profile in better agreement with predictions 
of galaxy formation models. This more fiexible parame- 
terization includes two parameters for the mean occupa- 
tion function of central galaxies and three for that of the 
satellite galaxies, which are described in detail below. 

The mean occupation function of central galaxies is a 
step-like function parameterized by 



(7Vee„(Af)> = i 



1 + erf 



log M - log Mn 



where erf is the error function 
erf (a;) = — ; 



flog A/ 



dt. 



(2) 



(3) 



There are two free parameters: Mmin, the characteristic 
minimum mass of halos that can host central galaxies 
above the luminosity threshold, and criogM, the width 
of the cutoff profile. Note that Mmin here can be in- 
terpreted as the mass of such halos for which half of 
them host galaxies above the given luminosity thresh- 
old, i.e., {Nccn{Mjnin)) = 0.5. It is n ot identical to thos e 
used in three-parameter models (e.g.. lZehavi et aLll2005D : 
(^ccn(-Mmin)) = 1/e for an exponential cutoff profile and 
{Nccn) changes from to 1 at Mmin for a sharp cutoff 
profile. Nevertheless, in all cases, Mmin characterizes the 
minimum mass scales of the host halos. In the three- 
parameter model, (A^con(M)), which uses one parameter 
Mmin, can also have a soft cutoff profile, e.g., an expo- 
nential profile as in IZehavi et al.l (|2005( ) . However, the 



two parameters Mmin and aiog m of the cutoff profile in 
equation ([2]) have a more physical meaning. To see this, 
we briefiy review the derivation of the form of (A^con(M)). 
The motivation for such a form is the relation between 
the central galaxy luminosity and the host halo mass pre- 
dicted by galaxy formation models. Based o n predictions 
of SA and SPH galaxy formation models, IZheng et al.l 
()2005l ) show that the distribution of the central galaxy 
luminosity Lc at fixed halo mass M can be described by 
a log-normal distribution, and here we write it as 



P(logL,|M) = 



1 



2'Kaioe- 



■ exp 



(logL,-log(Le(M)))2 



2<L 



(4)" 

In halos of mass M, the mean occupation function of 
the central galaxies above a luminosity threshold imin 
is an integration of equation (g]) over logic, and the re- 
sult turns out to have the same form as equation ([2]) but 
with the argument of the erf function being [logimin — 
log(Lc(M))]/(\/2criogL). If the mass range of the cut- 
off profile is not large so that we can approximate the 
mean luminosity (Lc(M)) of central galaxies in halos of 
mass M as {Lc{M)) cx M^, it is straightforward to see 
the meaning of the parameters in equation ^ — Mmin 
is the mass of halos in which the mean luminosity of 
central galaxies is the luminosity threshold Lmin, and 
the width of the cutoff profile is related to the scatter 
of the central galaxy luminosity in halos of mass M as 
ciogM = V^ciogi/p. Therefore, by studying the HODs 
of galaxy samples with different luminosity thresholds, 
we would learn the distribution (mean and scatter) of 
central galaxy luminosity as a function of halo mass, from 
the cutoff profile and the mass scale of {Nccn{M)). We 
follow the above explanation when interpreting our mod- 
eling results. 

According to predictions fron i iV-body, SA , and 
SPH galaxy forrn ation models (jKravtsov et al.l 12004 
IZheng et al.l l2005f) . the mean occupation function of 
satellite galaxies (-/Vsat(M)) approximately follows a 
power law at the high halo mass end with the slope 
close to unity. At lower mass, (iVsat(M)) drops steeper 
than the power-law extrapolation and can be parameter- 
ized by [(M — Mo)/Af{]", where Mq is the mass scale of 
the drop, M{ characterizes the amplitude, and a is the 
asymptotic slope at high halo mass. Applying the same 
cutoff profile of the central galaxies, we assume the fol- 
lowing form for the mean occupation of satellite galaxies, 
for M > Mq: 



(A^sat(M)) = - 



1 + erf 



log M - log Mn 

Clog Af 



/M 



Mo 



(5) 

For modeling the two-point correlation function of 
galaxies, we also need to know the second moment of 
the occupation number in addition to the mean. Central 
alaxies sim ply follo w the nearest-integer d istribution 
Berlind fc Weinbe73 [200l IZheng et al.l 12005). and for 
satellite galaxies, we assume a Poisson d istribution that 
is consistent with theoret ical predictions (jKravtsov et al.l 
12004 IZheng et all I2005D . The spatial distribution of 
galaxies inside halos is assumed to be the same as the 
dark matter that follows the NFW profi le, a reasonable 
assumption on scales we model (e.g.. iNagai fc Kravtsovl 
l2005t iMaccio et al.ll200il . Overall, our basic HOD pa- 
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rameterization has a total of five parameters — two for 
(A^con(M)) (Mmin and criogM) and three for (Arsat(M)) 
(Mo, M{, and a). With respect to the sim ple three- 
parameter model used in IZehavi et al.l ()2005[ ). our five- 
parameter model essentially introduces one additional 
parameter to {Nccn{M)) and one to (NsntiM)) to char- 
acterize the cutoff profiles at low halo mass. 

For theoretical calculations of galaxy two-point cor- 
relation func t ions, we adopt the method proposed in 
Tinker et al] (|2005D . which improves that presented in 
Zhend ( 20041 ) by incorporating a more accurate treat- 
ment of the halo exclusion effect. The method is cal- 
ibrated and tested using mock catalogs and can reach 
an accuracy of 10% or better in calculating galaxy two- 
point correlation functions. More specifically, we adopt 
the "fig-matched" approximation (see their Appendix B) 
for a more efficient calculation. 

4. HOD MODELING RESULTS FOR DEEP2 AND SDSS 
GALAXIES 

We perform HOD modeling of the projected two-point 
correlation function Wp{rp) for each DEEP2 galaxy sam- 
ple. The jackknife covariance matrices estimated from 
the DEEP2 data are noisy due to the size of the vol- 
ume probed, and we therefore use only the diagonal ele- 
ments in measuring x^- In calculating for each sam- 
ple, in addition to Wp^Vp) we also include the galaxy 
number density, obt ained by integra ting the observed 
luminosity function (jWillmer et al.l [20061 and assigning 
a 10% fractional uncertainty to the result. We ap- 
ply a Markov Chain Monte Carlo (M CMC; see e.g., 
iGilks. Richardson, fc Sp icgclhaltei1 ll996f ) method to ex- 
plore the parameter space. For the purpose of com- 
parison, we also m odel the SDSS galaxy samples in 
IZehavi et all (|2005D using the same HOD parameteriza- 
tion and MCMC method. The full jackknife covariance 
matrices are used for the SDSS samples. We perform a 
test with the Mr < —20.5 SDSS sample by only includ- 
ing the diagonal elements of the covariance matrix in the 
fit and find that the best-fit is higher, which reflects 
the fact that the likely positive covariance among val- 
ues of the two-point correlation function in adjacent bins 
is neglected. However, the marginalized distributions of 
HOD parameters are found to be pretty similar to those 
with the full covariance matrix. The test indicates that 
there are not likely to be large systematic uncertainties 
in the inferred HOD parameters for the DEEP2 samples 
for which we use only the diagonal elements of the co- 
variance matrices in the fits. 

Best-fit HOD parameters for DEEP2 and SDSS sam- 
ples are listed in Table [TJ We also list two derived param- 
eters in the table: Mi (not M[), the mass scale of a halo 
that can on average host one satellite galaxy above the 
luminosity threshold, and bg, the large scale bias factor 
of galaxies. In this section, we briefiy describe the mod- 
eling results and fits to Wp{rp), focusing on the DEEP2 
galaxies. Detailed inspections of individual HOD param- 
eters and comparisons between occupation properties of 
DEEP2 and SDSS galaxies are presented in the next sec- 
tion. 

4.1. Results for DEEP2 Galaxies 

The modeling results for the four DEEP2 luminosity 
threshold samples are shown in Figure [TJ In each of 



the left panels, lines from model predictions are plot- 
ted together with the Wp{rp) data points. All fits are 
excellent. With five parameters (one of which is largely 
determined by the mean galaxy density once the other 
four are specified), this success may seem unsurprising. 
However, it is worth noting that the combination of the 
CDM halo distribution (fixed by our assumed cosmol- 
ogy) and an HOD model cannot fit arbitrary functions 
— this is a physically constrained model with vastly less 
freedom than, say, a 4th-order polynomial. Our HOD fits 
self-consistently produce the two-point correlation func- 
tion both on small scales and on large scales. On small 
scales, the two-point correlation function reflects how 
galaxies are spatially distributed inside halos. On large 
scales (e.g., ^ 10/i~"'^Mpc), the galaxy two-point corre- 
lation function is simply the matter correlation function 
multiplied by a constant (i.e., the square of the scale- 
independent galaxy bias factor) . Our HOD model nicely 
reproduces the large-scale shape of Wp(rp), which indi- 
cates that the matter power spectrum we use has the 
correct shape on these scales. 

The central solid line corresponds to the best-fit model 
and is bounded by two solid lines, which are the envelopes 
of predictions from models with Ax^ < 1. As we do 
not have the full error covariance matrix, these should 
only be interpreted as indicative of the associated 1-cr 
uncertainty. The best-fit line is decomposed into contri- 
butions from the one-halo term (intra-halo galaxy pairs; 
dotted lines) and two-halo term (inter-halo galaxy pairs; 
dashed lines), which dominate on small and large scales, 
respectively. The transition from the one-halo term to 
the two-halo term causes an inflection in Wp{rp), lead- 
ing to departures from a pure power law that are ob- 
serve d for both low redshif t galaxies (e.g.. IZehavi et al] 
[20041: [Hawkins et al.' '2003') and high-redshift galaxies 
(e.g., [Ouchi ct al. 2005; Lcc et al. 2006). If the transi- 
tion scale is defined as the scale where the two contribu- 
tions to Wp{rp) are equal, we find that it increases from 
~0.4/i~^Mpc to ~0.6/i~^Mpc as the luminosity threshold 
of galaxies increases, a manifestation that more luminous 
galaxies reside in more massive halos. 

The mean occupation functions for the four DEEP2 
samples are shown in the right panels of Figure [TJ The 
mean occupation function for all galaxies in the sample 
is decomposed into that for central galaxies and that for 
satellites. For each mean occupation function, we plot 
the envelope from models with Ax^ < 1. It can be clearly 
seen that the mean occupation function shifts to higher 
halo mass as the luminosity of galaxies increases. In gen- 
eral, a low-mass cutoff in the satellite mean occupation 
function (which is usually neglected in the simple three- 
parameter model) is required by the DEEP2 Wp{rp) data 
for a good model fit. The best-fit models for the two 
brighter samples have a shallow cutoff profile in (Nccn)- 
For the two fainter samples, the inferred cutoff profile 
in (Nccn) can be either sharp or shallow, which leads to 
the kink seen in the envelope of (Nccn) for models with 
Ax^ < 1. The general trend that the central galaxy 
cutoff profile is better constrained for brighter galaxies 
reflects the fact that the halo mass function and bias fac- 
tor change steeply toward the high mass end, although 
the residual redshift-space distortion [because of Wp{rp) 
being derived from a finite projection] could artificially 
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TABLE 1 

Best-fit HOD Parameters for DEEP2 and SDSS Galaxy Samples'* 



Sample log Mmin 



flog M 



log Afo log AI[ 



logA/ii^ 



"9 



DEEP2 



Mb < -19.0 
Mb < -19.5 
Mb < -20.0 
Mb < -20.5 



11 64+°°* 

11 SQ+0.08 

12.07tH^9 
12.63l°- 



0.31 
0.30 
0.37 
0.82 



+0.19 
0.19 
+0.18 
0.19 
+0.17 
0.19 
+0.09 
0.09 



12 02"^"-"'^ 

il.OO_g gg 

Q no + l.SS 

^•'-'^-0.96 
o c:o+0.97 
».O»_0.94 



12 

J-^.O'_0.47 
13 02+° " 
1Q 97+0.06 

10.Z( _g Qg 

"+0.06 



13.56^ 



0.06 



0.89 
0.97 
1.08 
1.27 



0.23 
0.34 
+0.11 
0.11 
+0.05 
0.06 
+0.12 
0.12 



13 00+°-^° 
13 06+° °* 
i3.2«_(, (jg 

in co+o.oe 

1J.OO_0.06 



1.22 
1.32 
1.44 
1.47 



+0.01 
0.01 
+0.01 
0.01 
+0.01 
0.01 
+0.02 
0.02 



18 
18 
18 



8.4 
4.9 
7.9 



18 14.8 



SDSS 



Mr < 

Mr < 

Mr < 

Mr < 

Mr < 

Mr < 

Mr < 

Mr < 

Mr < 



-18.0 
-18.5 
-19.0 
-19.5 
-20.0 
-20.5 
-21.0 
-21.5 
-22.0 



11.35 
11.46 
11.60 



+0.07 
0.07 
0.06 
0.06 
+0.06 
0.06 

11 75+0.06 
^J^- "^-0.06 

12 02+°-°'' 
-0.05 
-0.05 
-0.10 
-0.10 



12. 3o: 

12.79; 
13.38j:°:i5 



14.22 



+0.13 



0.25 
0.24 
0.26 
0.28 
0.26 
0.21 
0.39 
0.51 
0.77 



+0.18 
0.16 
+0.18 
0.16 
+0.17 
0.17 
+0.16 
0.18 
+0.16 
0.16 
+0.16 
0.14 
+0.20 
0.22 
+0.22 
0.27 
+0.12 
0.18 



11.20 



+0.66 



-0.54 
.1+0.86 
-0.82 
3+0.78 
-0.96 
.1+0.68 
^-0.75 
oo + l.Ol 
11.3»_o.93 

11.84+°« 



10.59: 

11.49:* 
11.1 



12.40t°;l^ 
12.68t°:°° 

,+0.09 



0.42 
+ 1.15 
1.27 
+0.48 
0.51 

I4.00t°;i 



11.92 
13.94 



12.83 
13.01 
13.31 
13.58 
13.94 
13.91 
14.69 



0.08 
+0.09 
0.09 
+0.09 
0.09 
+0.06 
0.06 
+0.08 
0.08 
+0.56 
0.84 
+0.21 
0.18 



n OO+0.09 

-0.05 

1 02+°°'^ 
^■^^-0.05 
..+0.06 
-0.06 
..+0.08 
^-0.08 
1.12+°°'' 



1.06Z 
1.06^ 



1.15 
1.04 
0.87 



0.04 
+0.09 
0.10 
+0.55 
0.59 
+0.35 
0.52 



12.47 
12.70 
12.88 
13.06 
13.34 
13.59 
13.98 
14.47 
14.88 



+0.10 
0.10 
+0.08 
0.08 
+0.08 
0.07 
+0.08 
0.06 
+0.07 
0.07 
+0.06 
0.06 
+0.06 
0.06 
+0.07 
0.07 
+0.12 
0.12 



0.91 
0.95 
1.01 
1.03 
1.05 
1.15 
1.28 
1.52 
1.91 



+0.02 
0.02 
+0.02 
0.02 
+0.01 
0.01 
+0.01 
0.01 
+0.02 
0.02 
+0.01 
0.01 
+0.02 
0.02 
+0.05 
0.05 
+0.09 
0.08 



12 
12 
12 
12 
12 
12 
12 
11 
7 



12.7 
11.8 
8.6 
3.9 
6.2 
2.8 
2.3 
5.1 
1.0 



See § 3 for the HOD parameterization. Mass is in unit of h Mq. Error bars are for 1-a, which are derived from 
the marginalized distribution.'' These are derived parameters: Mi, the mass scale of a halo that can on average host 
one satellite galaxy above the luminosity threshold; bg, the large-scale galaxy bias factor.'^ This is the total number 
of data points (values of Wp plus the number density) used in the fitting. 



increase the shallowness of the cutoff profile (J. L. Tin- 
ker, private communication). 

We note that the constraints on the HOD models can 
be further strengthened by enforcing the condition that 
at any halo mass the mean occupation number for cen- 
tral or satellite galaxies in a sample with lower luminos- 
ity threshold is always higher than that for galaxies in a 
sample with higher luminosity threshold. A full imple- 
mentation of this would require modeling the data from 
all galaxy samples simultaneously, and there would be 
strong correlations among mean occupation functions of 
different samples. Here we choose to model each sample 
individually. To roughly account for the above condi- 
tion, we drop some models that are in apparent conflict. 
Specifically, the mean occupation functions of satellites 
for the Mb<— 19.0 sample from a small fraction of mod- 
els have a cutoff at much higher mass than those for the 
Ms <— 19.5 sample, which is not realistic. We therefore 
only keep models for the Mb<— 19.0 sample with a satel- 
lite cutoff mass of Mo<10^^h~^MQ for the following dis- 
cussions. We note that the above consistency condition 
is automatically satisfied by adopting the CLF method, 
with the HOD of each sample obtained by integrating the 
CLF over halo mass, although to reach a similar level of 
flexibility in the HOD we use here, more parameters need 
to be introduced in the CLF. 

The high-mass end slopes a for mean occupation func- 
;,iv,c. oofoiiifQo CT-Q fr^ninri i-n Kq n so+0.23^ 0.97^0}}, 



-0.34 



tions of satellites are found to be 

LOS^Qpg, and 1.27^0 12 in order of increasing luminos- 
ity threshold, where each value is quoted as the mean 
with error bars from the central 68.3% of the marginal- 
ized distribution (i.e., 1-cr range). Although the value 
for the most luminous sample is slightly larger, all val- 
ues are close to unity, which ag rees with predictions 
of galaxy formati on models (e.g.. iKravtsov et "all 120041 : 
iZheng et al.l[200l . 



We note that the best-fit HOD model for the bright- 
est DEEP2 sample {Mb < -20.5) has a larger than 
those of the other samples. Much of the difference be- 
tween the model and the data happens on small scales 
{rp < 0.5/i~^Mpc). On these scales, the raw measure- 
ments of Wp{rp) are affected by the undersampling of 
galaxies caused by the slitmask target selection algo- 
rithm, and they are correcte d by using the mo ck galaxy 
catalogs of lYan et al] (|2004l ) (jCoil et al.ll2006af ). In gen- 
eral, the correction increases the amplitude and the slope 
of Wp{rp) on small scales by a small amount. As a test of 
the robustness of our HOD results to this correction, we 
perform an MCMC run with the uncorrected •Wp{rp) data 
for this sample and find x^/d.o.f. — 0.75. The inferred 
HOD therefore seems to be robust against such a cor- 
rection — there is almost no change, except for a slight 
decrease in the number of satellite galaxies in low-mass 
halos to account for the slightly reduced amplitude and 
slope of the uncorrected Wp[rp). The slitmask correction 
could in principle be improved by using catalogs gener- 
ated according to the best-fit HOD model and iterating 
the correction and model fitting a couple of times. How- 
ever, given the robustness of the best-fit HOD under the 
adopted parameterization, we do not intend to perform 
such iterations in this paper. 

We also tried to fit the data with a more fiexible HOD 
parameterization, which uses a spline line for the satel- 
lite mean occupation function (similar to what is used in 
Fig. 19c of , Zehavi et al. ,20051 . We found that the HOD 
is adjusted to fit almost every feature in Wp{rp), and the 
modeling may result in solutions that are not seen in 
galaxy formation models. For instance, for the brightest 
Mb < —20.5 sample, the best-fit satellite mean occupa- 
tion function has an inflection and flattens out toward 
the low halo mass end, which allows it to fit the high 
amplitude and steep slope of Wpirp) on small scales that 
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Fig. 1. — HOD fits to tOp(rp) (left) and the mean occupation functions (right) for DEEP2 galaxy samples. In each Wp(rp) panel, 
overplotted on the data points are the range of lines predicted from models with Ax^ < 1 (the reduced from the best fit is labeled in 
each panel). The one-halo (dotted line) and two-halo (dashed line) terms from the best-fit model are also plotted for illustration. In each 
(N{M)) panel, the envelopes of mean occupation functions from models with Ax^ < 1 are plotted and the total mean occupation function 
(solid line) is decomposed into contributions from central (dashed line) and satellite (dotted line) galaxies. 



are vulnerable to slitmask corrections. Therefore, we 
conclude that this more flexible HOD parameterization 
might be "over fitting" the data, and that to robustly in- 
terpret the present data does not require an HOD form 
more flexible than our flve-parameter model. 

4.2. Results for SDSS Galaxies 

The modeling results for SDSS galaxies using the five- 
parameter model described in §3.2 are shown in Fig- 
ure [21 To be concise, only the results of the samples that 
are most relevant to the evolution connection between 
DEEP2 and SDSS galaxies (see § [B]) are plotted in this 
figure, although we perform HOD modeling for all SDSS 
samples with luminosity cuts ranging from — —18.0 
to Mr = -22.0 (see Table 1). 



Unlike the DEEP2 samples, the cutoff profiles of 
{Nccn{M)) for the SDSS samples are loosely constrained. 
This can be partly explained by noticing that in the mass 
range of 1011-101^/i~1Mq, both the halo mass function 
and the halo bias factor at z^l are steeper than those at 
z~0, and therefore the galaxy number density and the 
amplitude of the two-point correlation function at large 
scales lead to a better constraint in the cutoff profile for 
DEEP2 galaxies. We also note that complete freedom 
in CTiogM [eq. (2)] leads to a strong correlation between 
Mmin and txiogM in the sense that larger Mmin corre- 
sponds to larger ciog a/, which results in poor constraints 
on A/min. We assign a prior uiogA/ < 0.5, being conser- 
vativ e according to theoretical predictions (jZheng et al.l 
l2005f ). for SDSS galaxy samples with luminosity thresh- 
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Fig. 2. — Similar to Fig.[T] but for SDSS galaxy samples. 



old fainter than (A/;=-20.44). 

The transition scales from one-halo term to two-halo 
term for SDSS samples are around l/i~^Mpc, larger than 
those for DEEP2 galaxies. The cause of such a shift 
is that the one-halo term for SDSS galaxies appears to 
be shallower than that for DEEP2 galaxies, which again 
reflects the change in the slope of the halo mass func- 
tion in the mass range probed by DEEP2 and SDSS 
galaxies (note that, for the halo definition we adopt, 
at a given mass the comoving sizes of halos at the two 
epochs are the same). This also explains why the depar- 
ture of the two-point correlation function from a power 
law become s more prominent at highe r redshift (e.g. , 
2004; .Zheng. .2004 : .Kravtsov et all l2004t 



Zehavi et a 



Ouchi et al 



120051) ■ 



COMPARISONS OF THE HOD MODELING RESULTS 
FOR DEEP2 AND SDSS GALAXIES 




IQll 1012 1013 IQM 

M (h-iMg) 



In this section, on the basis of the MCMC results pre- 
sented in the previous section, we compare the HODs of 
DEEP2 and SDSS galaxies in detail. 

5.1. The Central Galaxy Luminosity as a Function of 
Halo Mass 



As mentioned in § 13.21 the low-mass cutoff profile of 
the mean occupation function of central galaxies encodes 
information about the distribution of central galaxy lu- 
minosities in halos with mass close to the cutoff mass. By 
modeling the galaxy two-point correlation functions for 
samples with different luminosity thresholds, we can de- 
rive the distribution of the central galaxy luminosity as a 
function of halo mass. If the central galaxy luminosity Lc 
in halos of fixed mass follows a log-normal distribution, 
as suggested by SA and SPH galaxy formation models 
(jZheng et al.ll2005l ). then the mass scale Mmin in equa- 
tion ([2]) is the mass of halos for which the mean central 
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Fig. 3. — Distribution of central galaxy luminosity as a function of halo mass. Left: Mean luminosity of central galaxies as a function of 
halo mass. Filled squares and circles are results for SDSS and DEEP2 galaxies, respectively. The error bars i n halo mass indicate the I-ct 
range of the distribution from MCMC runs. Dotted lines are calculated using the fitting formula proposed by Vale & Ostrikcr ( 2006) with 
Lo set to 2.8 X lO^ft-^L© (SDSS) and 4.3 X lO^fe-^L© (DEEP2) (see text for details). Right: Width of the distribution of central galaxy 
luminosity as a function of halo mass (see text for details) . 



galaxy luminosity (Lc) is simply the luminosity thresh- 
old of the galaxy sample (as shown in § 13. 2p . The left 
panel of Figure [3] shows the inferred {Lc)-M relation for 
DEEP2 galaxies and for SDSS galaxies. In the plot, ab- 
solute magnitudes are converted to luminosity in units 
of solar luminosity by using the Sun's abs olute magni- 
tudes 4.76 in r-band (jBlanton et al.ll2003bf l and 5.38 in 
i?-band^ (Johnson AB magnitude). 

The {Lc)-M relations look similar for the DEEP2 and 
SDSS galaxies except for an offset. At fixed halo mass, 
there is an offset of ^1.4 in (Lc) between the DEEP2 
and SDSS relations. However, given that the two sam- 
ples are selected in different rest-frame bands, this offset 
does not necessarily imply that z ~ 1 DEEP2 galaxies 
are more luminous than z ~ SDSS galaxies at a fixed 
halo mass. Alternatively, at fixed {Lc) there is an offset 
of ~ 1.6 in halo mass between the DEEP2 and SDSS 
samples. Similarly, the interpretation of such an offset 
is not immediately clear. In §[3] we attempt to establish 
a more meaningful evolutionary link between galaxies at 
the two epochs. 

Because of a large range in galaxy luminosity, the 
SDSS samples probe the {Lc)-M relation over 2 orders 
of magnitude in halo mass, from ~2 x 10^^h~^MQ to 
~3 X 10^^h~^MQ, as shown in the left panel of Figured 
We see that halos with higher mass host more luminous 
central galaxies at both z ^ 1 and z ^ 0. The plot shows 
that halos of Milky Way size (~2x IO^^/j-iMq) currently 
host central galaxies with mean lumino sity about 
(M: = -20.44, L* = 1.20 X IO^O/i'^Lq; iBlanton et aP 
I2003H ). Higher mass halos begin to host galaxy groups 
and clusters, so we expect that the {Lc)-M relation 
changes slope with halo mass. Indeed, the mean lu- 
minosity of central galaxies increases steeply with halo 
mass at the low-mass end and increases more slowly at 
the high-mass end, with the slope changing continuously 
from ~2.5 to ~0.3 over the mass range. 

The (Lc)-M relation from the SDSS modeling results, 

^ http:/ /www. ucolick.org/~cnaw/sun. html 



whose slope becomes shallower toward higher halo mass, 
is in general agreement wi th results f o r loca l galaxies 
derived by other methods. lYang et all (|2005| ) infer the 
(Lc)-M relation in bj band through associating halos 
with galaxy groups identified in the 2dFGRS, and they 
find (Lc) oc for halos of M < lO^^/i^iMpc and 

(Lc) oc Afi/^ for more massive halos. Using the X- 
ray masses of clusters/groups and the K-hand lumi- 
nosity of the brightest cluster galaxies in the Two Mi- 
cron All Sky Survey f2MASS). iLin fc Mohd ([200l find 
(Lc) oc for Af > 2 X lO ^^^M^. By rnodehn g 

the SDSS galaxy lensing data in iMcKav et all (|2001h . 
ICoorav &: Milosavlievi^ (|2005D obtain the {Lc)-M re- 
lation for M < IO^^A/r, which_ can be a pproximated 
as (Lc) oc M"-''^. Vale fc Ostrlkeil (|2004l ) present an 
empirical model for the relation between galaxy lumi- 
nosity and halo mass by matching the galaxy luminos- 
ity function and halo /subhalo mass function (see also 
IVale fc Ostrikeil I2006D and find that the mass luminos- 
ity re lation can be we ll appr oximated by a double power 
law. IVale fc Ostriked ()2006( ) advocate a fitting formula 
of the following form: 



(ic) = Lo 



{M/M,Y 



(6) 



with Afe = 3.7 X lO^/i-^AfQ, a = 29.78, h = 29.5, and 
k = 0.0255. With these parameter values, Lq is the mean 
galaxy luminosity in halos of mass 3.46 x \Q'^^h~^MQ. 
The relation is steep at low mass and becomes {Lc) oc 
^^0.28 jjjgij niass. We find that this fitting formula 
provides a good description of our modeling results of the 
SDSS galaxies if io = 2.8 x WH-'^Lq (see the dotted 
line for SDSS galaxies in the left panel of Fig. [3]). 

The DEEP2 galaxy samples we model probe the (ic)- 
M relation at z ~ 1 over a smaller mass range, from 
~4 X 10"/i-iAfQ to ~4 X IQ^'^h-'^MQ. Central galaxies 
of luminosity L^iMJ , = -20.7, L*g = 2.56 x lO^°h~^'^LQ, 
IWillmer et al.|[2b06( ). which is slightly above the highest 
luminosity threshold we probe, tend to reside in halos 
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a few times more massive than halos at z '--^ 0. At 
Af < 2 X IQi^/j-i^/jg^ tijg shape of the {Lc)-M relation 
at z ~ 1 is consistent with that at z ~ 0. At higher halo 
masses, {Lc) appears to increase with M slightly more 
slowly at z ~ 1 than at z ~ 0. The dotted line passing 
through the DEEP2 measurements in the left panel of 
Figure [3] is from the fitting formula in equation ^ with 
Lo set to be 4.3 x IO^/i-^Lq. 

In addition to the mean luminosity of central galax- 
ies, with our modeling results, we can also study the 
width (TiogL of the distribution of central galaxy lumi- 
nosity at a fixed halo mass [see eq. ([4])] . The information 
is encoded in the width aiogAf of the cutoff profile of 
(-^cen(Af)) and the local slope p of the {Lc)-M line — 
ciogi = pcriogAi/V^ (see I3.2p . The results of the width 
of the central galaxy luminosity distribution are shown 
in the right panel of Figure [Sj where a spline line passing 
through the {Lc)~M points is used to infer the local slope 
p for the SDSS and DEEP2 samples, independently. The 
error bar on each point indicates the l-cr scatter in aiog l 
(from its marginalized distribution) at the halo mass. In 
general, the scatter is not well-constrained by two-point 
correlation functions for low-luminosity samples, as the 
halo mass function and halo bias factor are not steep at 
the low-mass end. Therefore, we can remark only on gen- 
eral trends in the scatter. We find that at a fixed halo 
mass, the scatter uiogL in the central galaxy luminosity 
for z ~ 1 DEEP2 galaxies is about 10%-20% higher than 
that for z ~ SDSS galaxies, except for the brightest 
DEEP2 sample. For samples at both redshifts, there is 
a trend that ciogL decreases with increasing halo mass, 
although the significance from the modeling results is 
marginal. In the mass range we can probe, criogL changes 
from ~ 0.3 to ~ 0.1 for both SDSS and DEEP2 galaxies. 
The value of uiogL approaches a constant of 0.11-0.12 
at masses above 2 x 10^^h~^MQ for both samples. 

The decrease of uiogi with increasing halo mass (ex- 
cept for the last SDSS data point), i f true, is con s istent 
with theoret ical expectations (e.g., I Zheng et al.l 120051 : 
ICroton et al. 2006) in the halo mass range probed here. 
The last SDSS data point seems to show an increase but 
with large uncertainty. Although theoretical models also 
predict that aiog l should increase with halo mass again 
at higher mass, it is at a mass much higher than probed 
here (e.g., M > 5 x lO^^/i-^M©). We note as well 
that the observed int rinsic scatter of the TuUy-Fisher 
(iTullv fc Fisheil I1977D relation (e.g., iGiovanelli eFaLl 
[1992D, ~ 0.20^.35 mag, implies criogL =0.08-0.14 if one 
assumes that the circular velocity is a one-to-one indica- 
tor of halo mass. There is also a trend in the observed 
intrinsic scatter decreasing as the circular velocity in- 
creases, similar to what we see in our results. We spec- 
ulate that the reason for the relatively large scatter in 
central galaxy luminosity for lower mass halos may be 
that the distribution of major star formation epochs is 
broad for these halos. 

5.2. Relation Between Halo Mass Scales of Central and 
Satellite Galaxies 

In ij lS.ll we determined the mass scales of central galax- 
ies as a function of galaxy luminosity. We now focus on 
the halo occupation of satellite galaxies and the rela- 
tion between the halo mass scales of central and satellite 



galaxies. 

In our HOD parameterization, the parameter Mmin 
can be interpreted as the characteristic minimum mass 
of the halo hosting a central galaxy above the luminos- 
ity threshold: (A^con) = 0.5 at M^in [see eq.([2])]. For 
satellite galaxies, the characteristic halo mass scale Mi 
is defined as the mass of a halo that on average can 
host one satellite galaxy above the luminosity threshold: 
(7Vsat(Mi)) = 1 [note that Mi is different from M{ in 
eq.®]. In Figure [5] we show Afmin and Mi as a function 
of galaxy luminosity for DEEP2 and SDSS galaxies. In 
the left panel luminosity is relative to L, at each red- 
shift, and in the right panel luminosity is in solar units. 
For central galaxies, the data points are essentially those 
in Figure [3] viewed at a different angle. As mentioned 
in § 15.11 the DEEP2 and SDSS galaxies are selected in 
different rest-frame bands, such that inferring evolution 
between the samples is complex. Comparing galaxies rel- 
ative to at each redshift can roughly account for the 
different selections and the general dimming of galaxies 
since z '--^ 1. The plot shows that the host halos of z ~ 1 
galaxies are more massive than those of z ~ L* 
galaxies, a trend of downsizing in terms of the host halo 
mass of "typical" galaxie s with time since z ~ 1. 

In iZehavi et ail ((20051) . the HOD modeling for SDSS 
galaxies reveals a relation Mi ~ 23Mynin- Because of 
the slightly different HOD parameterization, Mmin and 
M l in our analysis do not correspond exactly to those 
in IZehavi et al.l ()2005l ). Our modeling results of the 
SDSS galaxies, as shown in Figure [H can be expressed 
as Ml - ISMmin. I Zheng et al. (2005) present Mi-M^in 
relations in SA and SPH galaxy formation models for 
samples of galaxies above varying baryonic mass thresh- 
olds and find that the scaling factors are about 18 and 
14 for the two models, respectively. Their definitions of 
M^in and Mi are more consistent with the ones used 
here, and this may be the reason that our results here 
ar e in better a g reeme nt with the theoretical predictions 
of iZheng et all (l2005h . 

The Mi-Afinin scaling relation for z 1 DEEP2 galax- 
ies is similar to that for z ~ SDSS galaxies, as shown 
in Figure m The scaling factor is about 16, close to the 
value for SDSS galax ies. Dissipationless simulations in 
iKravtsov et al.l ()2004f ) show a weak trend of the scaling 
factor decreasing towards high redshift, which can be un- 
derstood as the satellite galaxies not having enough time 
to merge with the central galaxies at higher redshift. 

For both data sets, the scaling factor of the Mi-Mmin 
relation has a trend of becoming smaller at the high- 
luminosity end. High luminosity satellite galaxies reside 
in massive halos. At any given redshift, massive halos 
form late and the accreted satellite galaxies do not have 
enough time to merge with the central galaxies. This 
effectively reduces the mass of the halos that on average 
host one satellite galaxy. Therefore, the effect of merging 
timescale may be the primary reason that the scaling 
factor decreases for very luminous galaxies. 

From Figure [H we can also see that at a fixed halo 
mass central galaxies are much more luminous than typ- 
ical satellite galaxies for both DEEP2 and SDSS samples, 
with the difference being larger toward lower halo mass. 
This is a manifestation of a CLF bump cause d by the cen- 
tral g alaxy seen in galaxy formation models () Zheng et al.l 
[2005h . 
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Fig. 4. — Mass scales of halos hosting central galaxies and satellites as a function of galaxy luminosity. Left: Mass scales as a function 
of L/Lt. Filled symbols connected by solid lines are the relation between the characteristic minimum halo mass Af„jin ^'nd the luminosity 
threshold, and open symbols are that between Mi and the luminosity threshold, where A/i is the mass of a halo that on average hosts one 
satellite galaxy above the luminosity threshold. Dotted lines are obtained by multiplying the solid lines by a factor of 18 (SDSS) and 16 
(DEEP2), respectively. Right: Same as the left panel, but the luminosity is expressed in units of solar luminosity. 

Figure [5] shows the satellite fraction as a function of 
the luminosity threshold for DEEP2 and SDSS galaxy 
samples based on the MCMC runs, where the luminos- 
ity threshold is expressed in units of L*. SDSS galaxies 
have a larger satellite fraction than DEEP2 galaxies. For 
example, ~20% of i > SDSS galaxies are satellite 
galaxies, compared to ^--^10% of i > DEEP2 galax- 
ies. For both DEEP2 and SDSS galaxies, the satellite 
fraction tends to decrease as the luminosity threshold 
increases. For L > L*, the satellite fraction as a func- 
tion of luminosity seems to drop more steeply. This trend 
with galaxy luminosity reveals that luminous galaxies are 
much more likely to be central galaxies in lower mass ha- 
lo s than satellite galax ies in higher mass halos. 

IConrov et al.l ()2006l ) predict the evolution of the lumi- 
nosity dependence of galaxy clustering for the redshift 
range < z < 5 by monotonically relating galaxy lumi- 
nosities to the maximum circular velocity of halos and 
subhalos in A^-body simulations. Their Table 2 lists 
the relevant HOD parameters at each redshift for var- 
ious luminosity threshold samples, with different num- 
ber densities. Our inferred HODs at z ~ 1 and z ^ 
are in good agreement with their non-parametric results. 
The dotted lines in Figure [5] are satellite fr actions calcu- 
lated from interpolating the parameters in IConrov et al.l 
((200l . which are consistent with our results in the rel- 
evant luminosity range. We note that satellite frac- 
tions, obtained here through HOD modeling of the two- 
point correlation function, can in principle be compared 
to those inferred through relating galaxy groups iden- 
tified in the surveys (e.g., iGerke et a l. 2005; Coil et al 
|2006b for the DEEP2 survey and [Bcrlind et al. 200 ^ 
IWeinmann et ani2006l for the SDSS survey) to dark mat 
ter halos for a consistency test. 

On large scales, galaxy pairs have contributions from 
central galaxies paired with each other. Being lower in 
number, satellite galaxies do not contribute as much to 
the two-point galaxy correlation function on large scales 
as on small scales. However, applications that relate the 




L/L. 

Fig. 5. — Satellite fraction as a function of luminosity threshold 
for DEEP2 and SDSS galaxies. The luminosity threshold is ex- 
pressed in units of L* at either redshift. Symbols show the mean 
value of the satellite fraction, and the error bars indicate the 1— cr 
range of the marginalized distribution from MCMC runs. The dot- 
ted lines are calculated from interpolating the HOD parameters in 
IConrov et al . ( 2006 ) , which are inferred through relating halos and 
subhalos in Af-body simulations to galaxies. 



5.3. Satellite Fraction 

A galaxy at a given luminosity can be either a central 
galaxy in a relatively lower mass halo or a satellite galaxy 
in a higher ma ss halo, with t he former being much more 
probable (Zhe ng et al.l|2005[ ). Although at a fixed lumi- 
nosity satellite galaxies are not dominant in number, they 
play a significant role in shaping the two-point correla- 
tion function of galaxies on small scales, where the one- 
halo term dominates. On these scales, galaxy pairs are 
composed of central-satellite and satellite-satellite pairs. 
Therefore, the two-point correlation function can pro- 
vide strong constraints on the overall fraction of satellite 
galaxies. 
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large-scale galaxy bias factor (relative to dark matter) 
to the halo bias factor often neglect this small contri- 
bution from satellite galaxies, which can systematically 
affect the results on the inferred mass scale and number 
density of host halos. In the Appendix we quantify the 
systematic error resulting from the "one galaxy per halo" 
assumption using a simple HOD model. 

6. EVOLUTIONARY CONNECTIONS BETWEEN DEEP2 
GALAXIES AND SDSS GALAXIES 

Thus far, we have compared different aspects of the 
halo occupations of DEEP2 and SDSS galaxies, treating 
the results as "static" observations of the galaxy pop- 
ulation at two different redshifts. We now attempt to 
establish an evolutionary link between galaxies at the 
two epochs to study the evolution of the (statistically) 
"same" galaxies with time. By "same" we mean a galaxy 
and its most massive progenitor. 

HOD modeling converts observed galaxy clustering 
measurements to a relation between galaxy properties 
and their host dark matter halos. The formation and evo- 
lution of the halos themselves are dominated by gravity, 
which is well understood from theory and simulations. 
In principle, we know how the halo population at an ear- 
lier epoch evolves to that at a later epoch. Given the 
relationships between galaxies and halos at two different 
epochs, the galaxy populations at those epochs can then 
be linked through the growth of their halos. 

In what follows, we implement this idea to study the 
evolution of central galaxies in the DEEP2 and SDSS 
samples. We first establish a relationship between 2; ~ 
halo mass and progenitor halo mass at z '--^ 1 in § 16.11 
Using this relation and the HOD modeling results, in 
S 16. 21 we link the z central galaxies and their progen- 
itor central galaxies at z ~ 1 in terms of luminosity to 
measure the luminosity evolution of galaxies since z ^ 1. 
The stellar mass of a galaxy is perhaps a more funda- 
mental physical parameter than luminosity and is more 
straightforward to interpret and compare to theoretical 
models. Therefore, in ^ \6.3\ using the same halo mass re- 
lationship, we attempt to link central galaxies at the two 
redshifts in terms of stellar mass to study the growth of 
stellar mass as a function of host halo mass. From this we 
are able to draw tentative conclusions on galaxy evolu- 
tion and star formation as a function of halo mass during 
the past ~7 billion years. Finally, in ^ \QAl using the stel- 
lar mass derived in § 16. 3[ we study how star formation 
efficiency depends on halo mass at the two redshifts. 

6.1. Typical Growth of Halos from z ^1 to z ^0 

To use our HOD results presented here, what we want 
to know is the mass of the z ^ 1 progenitor halo for each 
of the z ^ host halos. A dark matter halo merger tree 
that traces the assembly history of halos serves this pur- 
pose well. We use the PINOCCHIO (PINpointing Orbit- 
Cr ossing Collapse d Hierarchical O bjects) code developed 
by [Monaco et a l.' (2002a; see also [Monaco et alll2002bt 
iTaffoni et al.l 12002) to predict the assembly history of 
dark matter halos. PINOCCHIO is based on Lagrangian 
perturbation theory and can generate synthetic catalogs 
of halos that include mass, position, velocity, merger his- 
tory, and angular momentum. For the halo assembly 
history, predictions of PINOCCHIO are much more ac- 
curate than those based on the extended Press-Schechter 



forma lism fEPS: iPress fc Schechte"^ 11974 iBond et all 
I199H iLacev fc Colelll993D . Compared to iV-body sim- 
ulations, PINOCCHIO can predict the average mass as- 
sem bly history of halos to an accuracy of 10% or bet- 
ter ( Li et al.ll2005| ) with several orders of magnitude less 
computational time. Furthermore, unlike TV-body simu- 
lations, which need expensive postprocessing to extract 
merger trees, PINOCCHIO directly outputs the merger 
history of each halo. We perform four realizations of 
PINOCCHIO in a lOOMpc^ box and three realizations 
in a SOMpc"^ box to probe assembly histories of halos 
with mass above and below IO^^Mq, respectively. For 
these runs, the grid size is set to be 128^. The mass dis- 
tribution of the z ~ 1 progenitor halos of the z ~ halos 
is obtained in a series of z ^ halo mass bins. The code 
does not include evolution of subhalos within parent ha- 
los, but it suffices for our purpose here as we focus only 
on central galaxies. 

The results of the typical halo growth from z ^ 1 to 
z ~ are presented in Figure[5J:, which shows the relation 
between M, the mass of z ^ halos, and Mp^o, the mass 
of their z ~ 1 progenitors. The solid line represents the 
mean mass of the progenitors, and the two dotted lines 
indicate the boundaries of the central 68.3% distribution 
of the progenitor mass as a function of the z ~ halo 
mass. On average, lower mass halos grow earlier, in that 
more of their final mass is assembled by z ~ 1. The 
results show that a typical z ~ halo with M ^ 3 x 
lO"/i-iM0 (M - lO^^/i-iMg) has about 70% (50%) of 
its final mass in place at z ~ 1. We also use a merging 
tree code based on the (less accurate) EPS formalism and 
find that it generally predicts late growth for halos and 
leads to ~15% lower masses of the z ^ 1 progenitors in 
the mass range we consider. 

6.2. Luminosity Evolution 

The question we attempt to answer is, for galaxies of 
a given luminosity (e.g., L*) at z ~ 0, what luminosity 
did their z ~ 1 progenitor galaxies have? We answer this 
question using the HOD results above and the Mp^o-M 
relation shown in Figure [B]:. 

The Mpio-M relation links halos at z ~1 to those at 
z ~0. Together with the (Lc)-Af relations at these two 
redshifts derived from DEEP2 and SDSS galaxies, which 
are plotted in Figure [6li and Figure [Bh, it enables a con- 
nection between DEEP2 and SDSS central galaxies. The 
dashed lines across the four panels of Figure [6] illustrate 
how such a connection is established. Figure[6l> shows the 
resultant connection between the r-band luminosity of 
z ~ central galaxies and the i3-band luminosity of their 
z ^ 1 progenitors. For example, the z ^ 1 progenitors 
of the z - {M* = -20.44; Lr = 1.20 x W^Ii-'^Lq) 
central galaxies on average have a i?-band luminosity of 
Mb = -20.0 {Lb = 1.34 x W^Ji-'^Lq), less than at 
z ~ 1. 

The above connection is between luminosities in two 
different rest-frame bands at the two redshifts. A more 
interesting comparison may be between luminosities in 
the same rest-frame band. The best way to achieve such 
a goal is to have galaxy samples observed in the same 
rest-frame band, which is not the case for samples an- 
alyzed here and is hard to achieve in general at differ- 
ent redshifts. However, we can obtain a rough connec- 
tion by making a transformation between luminosities in 
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Fig. 6. — Connection between DEEP2 galaxies and SDSS galaxies. Panels (a) and (d) show the mean luminosity of central galaxies as 
a function of halo mass at z ~0 (SDSS) and z ~1 (DEEP2). Panel (c) shows the relation between the mass M of the z ~0 halos and the 
mass Mpro of their z ~1 progenitors. The solid line is the mean relation, and the dotted lines mark the central 68.3% of the distribution. 
Panels (a), (c), and (d) lead to a luminosity connection between DEEP2 galaxies and SDSS galaxies, which is shown in panel (fe). The four 
dashed lines illustrate how this connection is established (see text for details). 

stars have a shorter lifetime. For a stellar population 
with 0.2-1 solar metallicity, which forms its stars in a 
single instantaneous burst at redshift z =1.5, 2, 2.5, or 
3, its rest- frame i?-band luminosity would decrease by 
about 1.8, 1.4, 1.2, or 1.0 mag, respectively, from z = 1 to 
2 = 0. Using these values as an estimate, s.z ^ \ progen- 
itor central galaxy would have Lb ~2.6-5.3x 10^/i~^Lo 
(higher luminosity for stars forming at higher redshift) 
if they were passively evolving to z ^ 0. This is fainter 
than the z ^ central galaxies [Lb ^ 8.0 x lO^/i^^L©), 
which implies that additional stars must have been added 
to the system since z ~ 1, through star formation and/or 
galaxy mergers. 

6.3. Stellar Mass Evolution 

The connection between luminosities of local central 
galaxies and their progenitors provides interesting in- 
formation on the evolution of the stellar component of 



the t wo bands. We use the publicly-available kcorrect 
code ()Blanton et al.|[2003al ) to estimate the median Mb 
values (where B is a Bessell AB magnitude) for SDSS 
galaxies in a series of narrow Mr bins, to facilitate com- 
parison with Mb samples in the DEEP2 data. Using the 
SDSS DR4 d ata and the same m agnitude and redshift 
thresholds as IZehavi et al.l (|2005D . we find that the me- 
dian Mb value for galaxies at a fixed Mr is Mb — Mr + 1 . 
Our transformation implies that z ~ L* galaxies have 
a median Lb ~ 8.0 x IO^H^^Lq. On average their 
z ^ 1 progenitors have Mb = —20.0 (Fig. Ef)), i.e.. 
Lb ^ 1-34 X lO^°/i~^L0. It is not entirely surprising that 
their progenitor galaxi es are more luminou s , as st ars fade 
as they age. Using the lBruzual fc CharlotI ()2003D model, 
wc find that the amount of passive luminosity evolution 
at a fixed time interval depends on the age of the stellar 
population, in the sense that younger populations fade 
more. This is related to the fact that luminous massive 
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Fig. 7. — Similar to Fig.[6l but in terms of mean stellar mass. 

galaxies from z ~ 1 to z ~ 0. However, a compar- 
ison in terms of stellar mass would be more informa- 
tive and straightforward to interpret, as stellar mass is 
a more fundamental physical parameter than luminos- 
ity. As with luminosities, the ideal way to make this 
comparison is to construct stellar-mass-selected galaxy 
samples and to perform HOD modeling of the observed 
clustering of these samples, which we reserve for a future 
study. With the luminosity-selected samples used here, 
we use a simple estimation of stellar mass derived from 
the galaxy luminosity and color. While the correlation 
between stellar mass and galaxy luminosity can have a 
fair amount of scatter, we attempt here to extract use- 
ful information from the mean relation. The following 
analysis is presented as a proof of concept as to how one 
could attempt to study evolution of galaxies in terms of 
stellar mass and as a function of halo mass, and it serves 
only as a useful first-order approximation for analyses 
of stellar-mass-selected samples. With these caveats in 
mind, we proceed as follows. 
For the DEEP2 galaxies, stellar masses are derived by 



iBundv et al.l (|2006D for the subset for which iiTg -band 
imaging exists, assuming a Ch abrier stellar initial mass 
function (IMF; (|Chabrieiil2003l )). An empirically-derived 
relation between rest-frame UBV colors, redshift, and 
stellar mass is then used to estimate masses for the rest 
of the DEEP2 sample (C. N. A. Willmer, private commu- 
nication). For the SDSS samples, we estimate the stellar 
mass from the g—r col or and the r -band luminosity using 
the relation given bv lBell et ahl (|2003f) . The dominant 
source of uncertainty in estimatin g stellar mass using th e 
above method is the IMF of stars ffiell fc de Jon3l2001h . 
The "diet" Salpeter IMF (jBell et al.ll2003D is used for the 
SDSS galaxies. The effect of different IMFs is largely an 
offset in the estimated stellar mass. For example, stel- 
lar mass with the diet Salpeter IMF (Chabrier IMF) is 
70% (50%) that with the Salpeter IMF. We multiply the 
DEEP2 stellar masses with the Chabrier IMF by a factor 
of 1.4 to convert them to a diet Salpeter IMF, so that 
stellar masses from both surveys can be compared. 

The mean stellar mass of galaxies and the scatter are 
calculated in a series of narrow luminosity bins. We can 
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Fig. 8. — Mean stellar mass in 2 central galaxies and that in 
their 2^1 progenitors (top) and their ratio (bottom) as a function 
of the 2 ~ host halo mass. 

then convert the relation between central galaxy luminos- 
ity and halo mass to that between stellar mass of central 
galaxies and halo mass. With the present-day central 
galaxies and their progenitors at z ~ 1 linked through 
the Mp^o-M relation, we recast Figures [5] in terms of 
stellar mass, which is shown in Figure [71 Here Figure [7{> 
shows the relation between the stellar mass of z ~ cen- 
tral galaxies and that of the progenitors at z ~ 1. In 
Figure [21 the (arithmetic) mean stellar masses in central 
galaxies at z ~ and those of their 2: ~ 1 progenitors 
as a function of the present-day halo mass are plotted, 
which allows an estimate of the average growth of central 
galaxy stellar mass from z^ltoz'-^^Oasa function of 
halo mass. Limited by the luminosity range of the data, 
the halo mass range that can be probed is slightly greater 
than 1 order of magnitude. The bottom panel of Figure[Sl 
shows the ratio of the mean stellar masses in the z ~ 1 
progenitors and the z '--^ central galaxies as a function 
of the present-day halo mass, which is expected to be in- 
sensitive to the choice of the IMF. At 5 x 10"/i~^Mq, on 
average, a central galaxy had ^^20% of its present stellar 
mass in place at z ~ 1. The ratio gradually increases to 
~33% around 2 x 10^'^h~^MQ and is roughly constant 
up to the highest halo mass we can probe. For galaxies 
that reside in 2 x lO^^/i~^M0 halos (i.e., z ~ L* galax- 
ies), the increase in stellar mass from z ~ 1 appears to 
be consistent with the luminosity comparison presented 
in § 16.21 if most stars in their z ^ 1 progenitors do not 
form at a very high redshift and the luminosity is roughly 
proportional to stellar mass. As shown in Figure [71:, at 
z ~ 1, progenitor halos have already reached more than 
50% of their total present-day mass. Therefore, in the 
mass range probed here, halos have most of their mass 
assembled by z ~ 1 but have most of their central galaxy 
stars assembled / formed more recently. The mass scale of 
l-2x 10^^/i~^Mq for present-day halos appears to be a 
transition scale, below which a relatively larger fraction 
of stars have been added to the central galaxies during 
the period from z ~ 1 to z '--^ 0. 

In general, there are two processes that can add stellar 
components to central galaxies. One is star formation, 



and the other is galaxy merging, either merging of satel- 
lites with central galaxies or merging of central galaxies 
in different halos. The star formation includes that which 
occurs in the central galaxies and the star formation in 
the satellite galaxies that eventually merge onto central 
galaxies. If we know the halo occupation as a function of 
stellar mass at z ~ 1 , we can evolve the subhalo and halo 
population at z ~ 1 to z ~ using high -resolution N- 
body simulations or an analytic code fe.g.. lZentner et aLl 
|2005[ ). assuming no star formation. Such a calculation 
would give the amount of stars acquired through merg- 
ers. The difference between the evolved stellar mass oc- 
cupation and the true stellar mass occupation in halos of 
a given mass (inferred from z data) would tell us the 
amount of stars formed during the period from z '--^ 1 to 
z ^ 0. While the samples we model are not well suited 
for such a sophisticated analysis, it is a goal we plan to 
pursue in future work. 

The results presented in Figure [51 nevertheless allow us 
to place interesting limits on the amount of stars added 
to the present-day central galaxies by star formation and 
merging, in an average sense. As seen from Figures[6land 
[51 for z ~ halos with M = 5 x lO"ft.-Uf0, their z - 1 
progenitors have assembled ^70% of the z ^ mass 
while the progenitor central galaxies have ~20% of the 
stellar mass in place. For low-mass halos, it is reasonable 
to assume that the stellar mass in the progenitor halo 
is dominated by those from central galaxies (see, e.g., 
iZheng et al.ll2005l ). The progenitor halo therefore assem- 
bles the remaining ^-^30% mass by z ~ through smooth 
mass accretion and/or merging with one or more smaller 
halos. If the fraction of stellar mass inside this additional 
mass is the same as in the progenitor halo (i.e., smaller 
halos are assumed to have the same star-formation effi- 
ciency as the progenitor of the 5 x 10^^h~^ Mq halo), then 
at most ~9% of the final stellar mass can be gained by 
this process. That is, the stellar mass in the progenitor 
galaxy and that in the smaller halos that could merge 
with it or the mass that is accreted can only amount 
to ^30% of the final stellar mass in the z ~ central 
galaxy. Therefore, for central galaxies in z ^ halos 
with M = 5 X 10"/i-^Mq, -70% of the stars should 
form between z ~ 1 and z ~ 0. 

Similar reasoning can be applied to higher mass ha- 
los, but one should keep in mind that the contribution 
of stellar mass in satellites increases with halo mass. 
For z ~ halos with the highest mass we can probe 
(Af < 10^^h~^MQ, which host galaxy groups), satellite 
galaxies contribute a substantial amount of the stellar 
mass in the halo. The z ~ 1 progenitor of an aver- 
age halo of this mass has assembled ~54% of the total 
mass, and the progenitor central galaxy has ^33% of its 
stellar mass assembled. Neglecting star formation, the 
progenitor central galaxy can merge with other central 
galaxies (with their total halo mass amounting to the re- 
maining ^46% of the z ^ halo mass), and the total 
contribution to the stellar mass in the final z ~ central 
galaxy from all the merged progenitor central galaxies 
would be at most 60%. Merging of satellites to the final 
central galaxy would make an additional contribution to 
the stellar mass. We crudely estimate this contribution 
as follows. In a halo, satellite galaxies are usually fainter 
than the central galaxy, and in galaxy groups, the bright- 
est satellites are likely to be ~1.6 mag (~25%) fainter 
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Fig. 9. — Illustration of contributions to the average growth 
(since 2 ~ 1) of stellar mass of central galaxies residing in 2 ^ 
halos, as a function of the z ^ halo mass. See text for more 
details of the estimation based on the HOD modeling of galaxy 
clustering for the DEEP2 and SDSS galaxy samples. The top solid 
line would move up to be the dotted line after an upward correction 
of 25% for the DEEP2 stellar mass estimation (see text). 

than the central galaxie s (i.e., the luminosity gap; see 
iMilosavlievic et al] |2006| and references therein). This 
means that even if all the brightest satellites in halos 
that are assembled into a z ~ 10^^/i~^Mq halo are 
able to merge onto the final central galaxy, their contri- 
bution to the stellar mass is only about 25% of that from 
merging of central galaxies. Thus, the total stellar mass 
from merging of central galaxies and satellites may likely 
account for ~ 60% + 60% x 25% = 75% of that in the 
final central galaxies. The remaining ~25% would then 
be the result of star formation between z ~ 1 and z ^ 0. 

Using this same line of reasoning, we illustrate in Fig- 
ure [HI the different contributions to the z ^ central 
galaxy stellar mass in the halo mass range we are able 
to probe. In this plot, a linear interpolation in logarith- 
mic halo mass, based on the above limits for the low- 
and high-mass ends, is assumed to estimate the contri- 
bution from satellite mergers. Based on these crude es- 
timations, we find an interesting result: on average, a 
large fraction of stars in z ~ central galaxies resid- 
ing in low-mass halos formed since z ^ 1 (e.g., ~70% 
for 5 X 10^^h~^MQ halos), while only a small fraction 
of stars formed for central galaxies in high-mass halos 
(e.g., ^^25% for < 10^^h~^ Mq halos) over the same pe- 
riod. This trend seems to be a manifestation of the 
so-called downsizing - a pattern in which the sites of 
active star formation shift from high-mass galaxies at 
early times to lower-mass systems at later epochs (e.g., 
ICowie et al.]ll996l : lJuneau et ani2005f l. If the trend con- 
tinues to higher halo mass, beyond the luminosity range 
we probe in this paper, there may be no substantial star 
formation occurring between z 1 and z ~ in very 
massive halos. 

Are these estimates reasonable? We can compare them 
to values obt ained from the glob al cosmic star formation 
history (e.g. jMadau et aT] ll996). probed through difFcr- 

ent techniqu es. Using a recent compilation of data in 

iFardal et all (|2006f) . we estimate that about 40% of aU 



the stars in the present universe formed between z 1 
and z = 0. Based on stellar rnass es timates from the 
COMBO-17 survev. lBorchet all (|2006D find that the to- 
tal stellar mass density of the universe has roughly dou- 
bled since z ~ 1. These numbers are well bracketed by 
our values. 

The mass dependence of star formation that we find 
here can also be compared with other results. We find 
that the stellar mass growth rate of central galaxies de- 
pends on the host halo mass. Since halo mass and 
central galaxy stellar mass are closely related, this im- 
plies that the growth rate depends on stellar mass itself. 
iJimenez et alJ ()2005l ) use spectroscopic modeling to infer 
the star formation histories of individual SDSS galaxies 
and derive the mean star formation history as a func- 
tion of galaxy stellar mass. By integrating the fraction 
of stellar mass in the last 6.8 billion years as a function 
of stellar mass (sum of the dotted and all the solid lines 
in the top panel of their Fig. 4), we find that for z ^ 
galaxies with stellar mass ~ 2 x lO^^M©, ~70% of their 
stars formed between z 0.8 and z ~ 0. This fraction 
drops to ~40% for galaxies with a z ^ stellar mass 
of ~ 2 X IO^^Mq. Using th e star formati o n hist ories of 
SDSS galaxies estimated bv iPanter et alJ (I2006D . w hich 
improves on the modeling in I Jimenez et al.l ~( 20051 ). the 
above numbers become ~60% and ~25%. Our estimates 
show reasonable qualitative agreement with these results. 
It is worth noting that ou r method differs significantly 
from that of iPanter et al.l (j2006), who use stellar popu- 
lation modeling of z ~ SDSS galaxies with no reference 
to z ^ 1 galaxies, and do not use the clustering of either 
population. Our present analysis uses the population 
synthesis model only to obtain stellar mass-to-light ra- 
tios (not evolutionary ages), and our conclusions about 
stellar mass growth result from the HOD modeling of 
clustering at the two different rcdshifts. The agreement 
between the results from the two methods is therefore 
impressive. 

Our results can also be compared with studies of 
the evolution of star formation rates (SFRs) as a func- 
tion of stellar mass. Analyzing the AEGIS (All Wave- 
length Extended Groth Strip In t ernatio nal Survey) data 
out to z ~ 1, iNoeske et all ()2007al ) find that star- 
forming galaxies form a distinct sequence with a power- 
law depe ndence o f the SFR on the stellar mass (see 
also Noeske et al.l L2007b). Converted to the "diet" 
Salpeter IMF, their result can be written as SFR — 
3.6M0 yr-i X (M,tar/lOi°M0)O-67 for M.tar in the range 
of 1.4-14xlOi°M0 and 0.2 < z < 0.7. If we parame- 
terize the total stellar mass growth rate from both star 
formation and mergers as a power-law, Mstar = /(^)-^s'tar 
(i.e., with a time-dependent normalization and time- 
independent power-law index) , we can then solve for the 
mean relation averaged over < z < 1 utilizing the 
stellar mass connection between these two epochs (i.e.. 
Fig. [7l>). For z central galaxies with stellar mass 
in the range of 0.7-5 x IO^^Mq, we find that the mean 
growth rate of their stellar mass during the past '--^7Gyr 
is Mstar = 2.9M0 yr-1 x (M,tar/lOi°M0)O-67. Compared 
with the SFR in the star formation phase mentioned 
above, it is interesting to notice that the general growth 
has the same dependence on the stellar mass but with 
only a ~20% lower normalization, which suggests that. 
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for these galaxies, star formation is the dominant mode 
of stellar mass growth and the star formation phase on 
average may occupy a substantial fraction of the 7 Gyr 
interval. Modeling galaxy clustering data in more red- 
shift slices between z ^ 1 and z ~ would lead to a 
better understanding of the implication of the general 
growth rate. 

The numbers derived here should be taken as a first- 
order estimate, as we have performed an approximate 
calculation as a proof of concept. There are two limit- 
ing factors in the current data. As mentioned earlier, 
the samples used here are defined using galaxy luminos- 
ity, and modeling stellar mass-selected samples would be 
more appropriate for measuring the evolution of the stel- 
lar mass growth rate as a function of halo mass. The 
other complication is the different rest-frame color se- 
lections of the galaxies, in particular the effect that the 
DEEP2 samples used are not entirely volume limited for 
red galaxies (see § 12. ip . This could lead to an underes- 
timate of the mean stellar mass, as red galaxies have a 
higher stellar mass than blue galaxies at a given lumi- 
nosity. Therefore, we may be overestimating the stellar 
mass evolution from z ^ 1 to z ^ 0. With smaller, 
volume-limited DEEP2 samples at a slightly lower red- 
shift, we estimate that the fraction of galaxies missing 
in each DEEP2 sample is ^--^10% and we also find that 
the mean stellar mass is likely to be underestimated by 
20%-30%. This means that ~25% more of the stellar 
mass could have been in place in the z ^ 1 progeni- 
tors than the numbers we quote here, which would move 
the line in the bottom panel of Figure [5] (and the low- 
est line in Fig. ^ up by ^--^25%. Consequently, the top 
solid line in Figure [9] would shift upward to the place 
of the dotted line. The star formation contribution to 
the stellar mass growth at the lower (higher) halo mass 
end we probe would change from ~70% (^25%) to ~65% 
(~5%). However, despite these uncertainties, the general 
trend of the contribution to the stellar mass growth as a 
function of host halo mass remains the same. Therefore, 
the halo mass-dependent evolution trend inferred from 
the data is a robust result. 

6.4. Evolution of the Star Formation Efficiency 

With the relation between stellar mass and halo mass 
at two epochs in hand, we now present the evolution of 
the star formation efficiency and discuss the implications. 

Assuming that the baryon fraction in halos equals the 
mean baryonic fraction fb in the universe, for which we 
adopt the value 15.7% (see § 13. ip . then stellar mass 
can be used to study the star formation efficiency as 
a function of halo mass at z and z ^ 1. Specifi- 
cally, we associate a halo of mass M with a baryon mass 
Mb = fbM and define Mstar/Mb as the star-formation 
efficiency, where Mgtar is the stellar mass of the central 
galaxy in the halo. This star formation efficiency is, in 
essence, the value integrated from the redshift of forma- 
tion to the epoch in consideration (z ~ 1 for DEEP2 
galaxies and z ~ for SDSS galaxies) and reflects the 
fraction of baryons associated with halos that have been 
converted to stars in the central galaxies by the given 
epoch. In practice, the baryon fraction in a halo may 
not equal the global fraction, due to various physical pro- 
cesses. For example, in a low-mass halo there may be a 
smaller fraction of baryons accreted because of the shal- 
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Fig. 10. — Mean fraction of baryons in host halos that are 
converted into stars in central galaxies for z ~ (SDSS) and z ~ 1 
(DEEP2) galaxies. 

low potential well. Therefore, more accurately, the halo 
star formation efficiency defined here is the true star for- 
mation efficiency times the baryon accretion efficiency of 
the halo. Nevertheless, it is still an instructive quantity 
that can be interpreted as an apparent star formation 
efficiency. We note that by definition this efficiency does 
not include stars in satellite galaxies, whose contribu- 
tion increases with halo mass (about a few percent at 
5 X IO^/i-^Mq and -40% at lO^^h-^MQ, according to 
the estimates in 16. Sp . 

Figure [TO] shows this star formation efficiency, reflect- 
ing the fraction of baryons that are converted into stars 
at the two epochs. Both lines are calculated from the ra- 
tio of the arithmetic mean stellar mass of central galax- 
ies in halos of a given mass to the average baryon mass 
associated with these halos. At z '--^ 0, there is a char- 
acteristic mass scale of halos, ~ 6 x lO^^ft.~^M0, where 
the average conversion efficiency from baryons to stars 
in central galaxies reaches a maximal value. Even at this 
peak, only ~27% of the baryons in the halo are converted 
into stars of the central galaxy. The star formation ef- 
ficiency drops steeply at lower halo masses and declines 
more slowly at higher halo masses. 

At z ~ 1, our inferred average conversion efficiency 
has a similar behavior to that at z ~ 0, but the overall 
conversion efficiency decreases and the peak is reached at 
a larger halo mass {^10^'^h~^MQ) with a value of ^--^12%. 
While the absolute value of the efficiency depends on 
the IMF, the trend as a function of mass and thus the 
c haracteristic ha l o mas s should not. 

IShankar et al.l (j2006[ ) infer the star formation efficiency 
at z ~ by matching the observed galaxy stellar mass 
function with the theoretical halo mass function. Their 
result (see the dashed line in their Fig. 5) is in a good 
agreement with ours on the trend of the efficiency with 
halo mass; our values are slig htly lower than the i rs but 
are well within the error bars. iMandelbaum et al] ()2006f ) 
perform HOD modeling of galaxy-galaxy lensing for stel- 
lar mass-selected SDSS galaxy samp les. T he stellar mass 
in their study assumes the Kroupal (|2001h IMF, which is 
about 30% lower than the diet Salpeter IMF we adopt 
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here. With the IMF difference corrected, our estimate 
of the dependence of the z ~ star formation effi- 
ciency on halo mass agrees well with their results (see 
their Table 3). After correcting for differences in the 
adopted IMF and halo definitions, in the relevant halo 
mass range our re sults at 2 ~ and z ~ 1 also agree with 
those reported bv lHevmans et all (|2006D based on a weak 
gravitational lensing study of the Hubble Space Telescope 
GEMS (Galaxy Evolution from Morphologies and SED) 
survey. It is worth emphasizing that our results follow 
from the mass distribution of CDM halos and the galaxy 
assignment required to reproduce the observed cluster- 
ing, while the galaxy lensing results are from direct mea- 
surements of the dark matter halos through their weak- 
lensing effects. The agreement between our results and 
the lensing results is therefore, again, impressive. 

What determines the halo mass scale where the aver- 
age conversion efficiency reaches a maximum? It could 
be sta r formation feedback or preheating. I Mo et ahl 
()2005[ ) argue that present day halos less massive than 
lO^^/i~^M0 were embedded in pancakes of M ~ 5 x 
lO^^/i~^M0 at z ~ 2, whose formation heats and com- 
presses the gas, leading to a cooling time longer than the 
age of the universe at z < 2. Therefore, in halos below 
IO^'^H^^Mq, there is not much cold gas available for star 
formation, which leads to the drop in the conversion ef- 
ficiency. Our z ^ mass scale is roughly consistent with 
their prediction. On the other hand, our results indi- 
cate that this halo mass scale shifts to a higher mass at 
z 1. If true, preheating alone may not be sufficient to 
explain this mass scale. The intensity of star formation 
and the subsequent feedback as a function of halo mass 
and redshift may also be an important factor. The shift 
of the peak of star formation efficiency from high-mass 
halos to low-mass halos with time can be regarded as 
another manifestation of the downsizing pattern seen in 
star-forming galaxies. For both DEEP2 and SDSS galax- 
ies, the drop of conversion efficiency above the character- 
istic halo mass could be in part due to the fact that we 
only consider stars in central galaxies, and in high-mass 
halos the stellar mass contributions from satellite galax- 
ies can be substantial. In addition, in high-mass halos 
gas accretion becomes less efficient because of the high 
virial temperature. We note again that the stellar mass 
results presented here are shown mainly as a proof of 
concept and any conclusions based on these should be 
regarded as tentative. 

7. SUMMARY AND FUTURE PROSPECTS 

We perform HOD modeling of the projected galaxy 
two-point correlation functions Wp{rp) for luminosity 
threshold samples of DEEP2 and SDSS galaxies at z 1 
and z ~ 0, respectively. The HOD modeling, which con- 
verts galaxy pair statistics to relations between galaxy 
properties and dark matter halos, reproduces well the 
galaxy correlation functions at these two redshifts, in- 
cluding the rise on small scales seen at z ^ 1. 

We infer the relationship between central galaxy lu- 
minosity and halo mass at z ~ and z ~ 1. We find 
that at both redshifts the mean central galaxy luminosity 
(Lc) increases with halo mass. In low-mass halos below 
~ 10^^ft.~^Mo, (Lc) increases more rapidly with halo 
mass, and the scatter in the central galaxy luminosity 
in halos of fixed mass appears larger, possibly reflecting 



the broad distribution of major star formation epochs in 
these halos. We find evolution in the HOD in that galax- 
ies at a given halo mass are ~1.4 times more luminous 
at z ~ 1 and for a given luminosity the halo mass is ~1.6 
times greater at z ^ 0, with the caution that luminosities 
of DEEP2 and SDSS are in different rest-frame bands. 
In addition, central galaxies are found to be in halos 
a few times more massive at z ~ 1 than at z ~ 0. 

For both DEEP2 and SDSS galaxies, there exists a 
scaling relation between Myain, the characteristic mini- 
mum mass of the halo that can host a central galaxy 
above a given luminosity threshold, and Mi, the mass 
scale of a halo that on average is able to host one satel- 
lite galaxy above the same luminosity threshold. There 
is little difference between the scaling relation at z ~ 1 
(Ml ~ 16Mmin) and that at z ~ (Mi ~ ISM^in). 
The fraction of galaxies that are satellites decreases with 
increasing galaxy luminosity at both redshifts. At a 
fixed luminosity threshold (in units of L*) SDSS galax- 
ies have a larger satellite fraction than DEEP2 galaxies, 
e.g., ~20% versus ~10% for L > galaxies. 

The ultimate goal in modeling galaxy clustering at dif- 
ferent redshifts is to learn about galaxy formation and 
evolution. These HOD modeling results provide us with 
relationships between galaxies and dark matter halos at 
2 ~ 1 and z 0^ spanning half the history of the uni- 
verse. We use the typical growth of dark matter halos as 
determined from simulations to establish an evolutionary 
link between DEEP2 and SDSS galaxies and extract in- 
formation about galaxy evolution over the last 7 billion 
years. We establish such a connection, on average, by 
linking z central galaxies to z ~ 1 central galaxies 
residing in the progenitor halos of the z ~ halos. We 
relate the luminosities of z '--^ central galaxies to those 
of their 2 ~ 1 progenitors; however, the interpretation 
is complicated somewhat by the different rest-frame se- 
lection bands of SDSS and DEEP2 galaxies. As a proof 
of concept, we use stellar masses determined by galaxy 
color and luminosity to roughly estimate the evolution of 
the dependence of stellar mass on host halo dark matter 
mass. At z 1 the mean fraction of baryons that have 
been converted into stars is below ~15%, peaking at a 
halo mass of IO^^H'^Mq. At z ~ this fraction is be- 
low ~30% and the peak of the star formation efficiency 
is at a lower halo mass of ~ 6 x 10^^h~^MQ. We find 
that, on average, the majority of stars in z central 
galaxies in low-mass halos (a few xl0^^h~^ Mq) formed 
between z ~ 1 and z ~ 0, while only a small fraction 
of stars formed in central galaxies of halos as massive as 
~ 10^^h~^MQ. This refiects the downsizing pattern seen 
for star-forming galaxies as a function of stellar mass, 
shown now as a function of their host halo mass. 

The stellar mass results presented here are prelimi- 
nary, as the galaxy samples we model are defined by 
luminosity, not stellar mass. In addition, the SDSS and 
DEEP2 surveys have different rest-frame selection func- 
tions, which complicates comparing the results at the 
two redshifts. Ideally we would like to have stellar mass- 
selected galaxy samples at different redshifts from sur- 
veys that have the same rest-frame selection functions. 
While our results apparently agree with other studies in 
many aspects, we regard our conclusions on the details of 
the stellar mass evolution as tentative. Nevertheless, the 
halo mass-dependent evolution trend wc infer appears to 
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be robust, and the exercise serves as a starting point and 
guide for future work. 

HOD modeling of galaxy clustering for stellar mass- 
selected samples at different redshifts would be a great 
tool for gaining insight on galaxy formation and evolu- 
tion, and we plan to pursue this in future work. Un- 
like the simple estimates presented in this paper, which 
are based on a mean relationship between the final halo 
mass and progenitor halo mass, a comprehensive pro- 
gram would be more sophisticated. Halos at high red- 
shift, populated with galaxies as a function of stellar 
mass in accordance with the HOD results, would serve 
as initial conditions. Then, assuming no star formation, 
these halos would be evolved forward in time to lower 
redshift, following the merging histories and satellite dy- 
namics determined from simulations or analytic models. 
This method is similar to the semi-analytic galaxy forma- 
tion model in merging of halos and dynamically evolving 
satellites, but here the initial conditions are set by the 
HOD modeling results. This method differs from semi- 
analytic models in that only dark matter evolution and 
dynamics, not baryon physics, are involved in the cal- 
culation. The difference between the evolved HOD and 
the HOD inferred from galaxy clustering at low redshift 
would then provide a wealth of information on star for- 
mation as a function of halo mass. For example, the 
contributions to stellar mass in central galaxies can be 
separated — what fractions are from star formation, from 
central galaxies in other halos that merge with them, and 
from minor mergers with satellite galaxies. We would 
also infer the average star formation history in satellite 
galaxies. If galaxy clustering is modeled at a series of red- 
shifts, then the average continuous star formation history 
and the average stellar mass growth history from other 
modes (i.e., major and minor mergers) as a function of 
halo mass can be determined. Such empirically-derived 
results would shed light on the physics of galaxy for- 
mation, and comparisons with predictions of galaxy for- 
mation models would provide stringent tests for galaxy 
formation theory. An exam ple of the effort starting along 
some of these Hues is Wh ite et al.l (|2007[ ). which uses 
iV-body simulations in combination with HOD modeling 
of clustering to put constraints on merging of luminous 
galaxies from z ^ 0.9 to z ^ 0.5. 

Additionally, the envisioned program may also help 
test the role of environment on galaxy formation and 
evolution. One of the main assumptions in the current 
HOD framework is that the galaxy content in halos de- 
pends only on the halo mass and is statisticall y indepen- 
dent o f the halo's larger scale environment. iGao et all 
(|2005[ ) find in simulations a signature of age-dependent 
halo clustering (the so-called halo as sembly bi as, see also 
Croton et all 120071; iJing et all 120071: iHarker eFaLl 1 2006 : 



Wechsler et al. ll2006HWetze" et al. 1 120071: IZhu et"alll2Ci06l : 
Gap fc White! [2007f ) . In general, the environmental ef- 
fect on halo clustering is essentially negligible for halos 
above the nonlinear mass scale and becomes stronger for 
lower mass halos. With the method outlined above we 
may be able to test the effect of large-scale environment 
on galaxy evolution. It is likely that the environment- 
dependent halo clustering does not play a significant role 
in interpreting the observed galaxy clustering at high red- 
shift due to the fast drop of the nonlinear mass (e.g., from 
z = 1 to z = 0, the nonlinear mass increases by a factor 



of ~25, while halos only grow by a factor of about 2 in a 
large range of mass). When we passively evolve a high- 
redshift HOD, which is insensitive to the environment, to 
lower redshift, the environmental effect during the evo- 
lution can be fully taken into account. Comparing then 
the evolved HOD (which accounts for the environmen- 
tal effect) with the HOD inferred from the low redshift 
data (which ignores the environmental effect) may pro- 
vide insight on how environment shapes the HOD and 
determines properties of galaxies in low-mass halos, in- 
creasing the constraining power of galaxy clustering on 
galaxy formation processes. 
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Fig. 11. — Left: Comparison between the minimum mass M^nin of the host halos under the assumption of one galaxy per halo (i.e., no 
satellite galaxies) and the true value at a measured galaxy bias factor for galaxy samples at different redshifts. The bottom panel shows 
the ratio of the Mmin from one galaxy per halo to the true value. Right: Same as the left panels, but for comparisons between the inferred 
number densities of halos (see text for details). 
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APPENDIX 

BIAS RELATED TO THE ASSUMPTION OF "ONE GALAXY PER HALO" 

Given a measured two-point correlation function for a galaxy population, HOD modeling is the ideal path to infer 
information about the host halos. However, this involves non-trivial modeling efforts, and therefore approximations 
are often used in the literature. One commonly-used method assumes that there is one galaxy per halo and connects 
galaxies to halos through the large-scale bias factor. This assumes that all galaxies in the sample are central galaxies, 
while in reality a small fraction are satellite galaxies. Using a simple HOD model, we quantify the resulting systematic 
bias in the inferred halo mass and the duty cycle of galaxies from this assumption. 

We consider galaxy samples defined by a lower threshold in luminosity or stellar mass. The usual procedure associated 
with the "one galaxy per halo" assumption is as follows. The large-scale bias factor of galaxies is obtained from the 
ratio of the measured galaxy two-point correlation function and the theoretical matter two-point correlation function 
under the assumed cosmology. This large-scale bias factor is the same as that o f halos above a mass threshold. Since 
for an assurned cos mology the latter is a known function of mass threshold (e.g.. Mo fc Whit3fl996l : ISheth et al.l[20nl 
[Tinker et al]|2005D . the halo mass threshold A/,nin can be determi ned. The number densi t y of halos above th e mass 
threshold Mmin can be inferred from the halo mass function (e.g.. ISheth fc TormenI [l99l iJenkins et al.ll2001[ ). The 
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ratio of the observed galaxy number density to the inferred halo number density is used to determine the duty cycle 
of galaxies. More exactly, the "one galaxy per halo" assumption should be called "/ galaxy per halo," where / is a 
positive constant smaller than unity representing the duty cycle. 

If satellite galaxies are taken into account, using the same galaxy bias factor measured from the data, a different 
A^min would be inferred. To see how large the difference is, we assume a simple mean occupation function to represent 
the "truth;" this model is the sum of a step function for central galaxies and a power law for satellite ga l axies , 
(N(M)) = 1 + M /Mi for M > M^in- Following galaxy for mation model pred ictions (e.g.. iKravtsov et all 12004 
IZheng et al.ll2005l ) and HOD modehng results (e.g., this paper; IZehavi et al]|2005[) . we use Ml = 20iVfmin. The above 
(N) can have an overall normalization factor smaller than unity (duty cycle factor), and none of our results below 
change. 

The bias factor b for a given galaxy sample (in which all galaxies reside in halos of mass M > Mj^in) can be calculated 

as 

r+°° fin I /■+°° dn 

6(>AW-y^^ dM—{N{M)}b,{M)^ j^^ dM—{N{M)), (Al) 

where dn/dM is the halo mass function and bh{M) is the bias factor of halos of mass M (following iSheth &: TormenI 
I1999D . Solid lines in the top left panel of Figure [TT] show the relation between Mmin and b{> Mmin) at different 
redshifts. Assuming one galaxy per halo, the bias is calculated by replacing (iV(M)) with unity in equation (|A1[) . 
Dotted lines in the top left panel of Figure [TT] show this relation. At a given bias factor (measured from the galaxy 
two-point correlation function), the value of M„iin with the assumption of no satellite galaxies {dotted lines) is always 
higher than the "true" value {solid lines) because massive halos (with high halo bias) are weighted more as the number 
of satellites increases with halo mass. The bottom left panel of Figure [TT] shows the ratio of the Mmin from the "one 
galaxy per halo" method to the true value at each redshift, as a function of the galaxy bias. When the bias is high, the 
corresponding halos are in the exponential tail of the mass function and the weight from the satellite galaxies residing 
in massive halos is less significant, therefore, the value of Mmin inferred from the approximate method approaches its 
true value. However, the difference increases for lower galaxy biases, corresponding to the shallower part of the halo 
mass function, and at higher redshift. For instance, at z = 0, for a galaxy sample with b = 1.2, the estimated A/min 
can be a factor of 2 higher than the true value. 

The difference in the inferred values of Mmin leads to a difference in the inferred number densities of halos above 
Mmin- The halo number density nh{> Mmin) is calculated as 

nh{> M„un) = / dM-^. (A2) 

The right panels of Figure [Tl] compare the inferred (from "/ galaxy per halo" ) and the true number densities of halos 
as a function of galaxy bias. For a given galaxy bias, the halo number density derived from the approximate M^i^ is 
always lower than the true value because without satellites in massive halos, a larger value of M^i^ is needed to reach 
the given bias factor. As with Mmin, the estimated halo number density approaches the true value for highly biased 
galaxy samples, while for low bias factors, the difference can again be substantial. For example, at z = 0, the number 
density of host halos for galaxies with b = 1.2 is underestimated by a factor of 2. As a consequence, the duty cycle 
(the fraction of halos hosting galaxies of that type) is overestimated by using the underestimated halo number density. 

For galaxy samples defined by a lower threshold (e.g. in luminosity). Figure [TT] can be used to better estimate the 
minimum host halo mass and number density from the measured galaxy bias value, instead of assuming one galaxy 
per halo. 
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